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The motions of a passive scalar a in a general high-frequency oscillating flow are studied. 
Our aim is threefold: (i) to obtain different classes of general solutions; (ii) to identify, 
classify, and develop related asymptotic procedures; and (iii) to study the notion of drift 
motion and the limits of its applicability. The used mathematical approach combines a 
version of the two-timing method, the Eulerian averaging procedure, and several novel 
elements. Our main results are: (i) the scaling procedure produces two independent di- 
mcnsionlcss scaling parameters: inverse frequency 1/lu and displacement amplitude S; (ii) 
we propose the inspection procedure that allows to find the natural functional forms of 
asymptotic solutions for 1/uj — > 0, 5 — > and leads to the key notions of critical, sub- 
critical, and supercritical asymptotic families of solutions; (iii) we solve the asymptotic 
problems for an arbitrary given oscillating flow and any initial data for a; (iv) these 
solutions show that there are at least three different drift velocities which correspond to 
different asymptotic paths on the plane (l/u;,<5); each velocity has dimensionless mag- 
nitude O(l); (v) the obtained solutions also show that the averaged motion of a scalar 
represents a pure drift for the zeroth and first approximations and a drift combined with 
pseudo- diffusion for the second approximation; (vi) we have shown how the changing of 
a time-scale produces new classes of solutions; (vii) we develop the two-timing theories 
of a drift based on both the GiM-theory and the dynamical systems approach; (viii) 
examples illustrating different options of drifts and pseudo-diffusion arc presented. 



1. Introduction 

The transport of a scalar field a by an oscillating velocity field u represents a classical 
problem of fluid dynamics. One can find a number of models, physical effects, and applica- 
tions in|Taylor (1953)||Moffatt (1983)||Hydon k Pedley (1993)||Frisch (1995)l|Moffatt (1998)] 



Magar k Pedley (2005) Leal (2007) The transport of a scalar represents a paradigm 



the studies of turbulence (see Sreenivasan (1991) Warhaft (2000) Balk (2002)). It is 



problem for the use of the dynamical systems approach in fluid dynamics (see Aref (1984) 
Ottino (1989)||Solomon, Lee, and Fogelman (200T)]|Samelson k Wiggins (2006)[) a nd for 



well-known that even in purely oscillating flows the motion of a material particle con- 
sists of two parts: oscillating and non-oscillating. The latter is known as a drift, see 



Stokes (1847) 


Maxwell (1870) Lamb (1932) Longuet-Higgins (1953) Darwin (1953) 


Hunt (1964) 


Batchelor (1967)| 


Schlichting (1979)[ Lighthill (1978a) Andrews k Mclntyre (1978)[ Craik (1982)[ 


Grimshaw (1984) 


Craik (1985) Benjamin (1986)| Eames k Mclntyre (1999) Barnes, Belcher, and Hunt (1994) 


Different notions of a drift are used in physics, see e.g. Vergassola k Avellaneda (1997) 



Chierchia k Gallavotti (1994) 



In our paper we consider a high-frequency asymptotic problem of a scalar trans- 
port by a given oscillatory velocity field. A characteristic length L for both u and 
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a is arbitrary. We use a simple and straightforward formulation of the problem that 
allows to study the general properties of solutions without the consideration of par- 
ticular applications. The evolution of a is exploited for the study of the concept of a 
drift motion and the limits of its applicability; this approach allows a great degree of 
flexibility and has not been used before. We employ the mathematical method which 



comprises the two-timing approach (Nayfeh (1973)1, the fast-time averaging operation, 



and several novel elements; we call it two-timing-and-averaging method (TTAM). The 
interplay of the two-timing form of solutions and a fast-time-averaging operation was 
exploited by many authors. The well-known applications in classical mechanics belong 
Kapitza (1951a)[ |Kapitza (1951b)| |Landau fe Lifshits (2000)[ A simplified version of 



to 



the two-timing method has been extensively used in mechanical engineering and is 
known as vibrational mechanics, see e.g. Blekhman (2000)| |Blekhman (2004) However 



Kapitza (1951a)||Kapitza (1951b)||Landau fc Lifshits (2000)[|Blekhman (2000j[|Blekhman (2004) 



did not use (and even deliberately denied) the underlying asymptotic nature of the 



method. This nature has been emphasised and exploited by Vladimirov (2005) Yudovich (2006) 



Vladimirov (2008) in a particular form of TTAM used in this paper. 
Let us outline here the content and the main results: 

In Sect. 2 we formulate the problem, introduce two-timing and dimensionless versions 
of the governing PDE. The two independent dimensionless scaling parameters are inverse 
frequency 1/cj and displacement amplitude 5 of oscillations. Our asymptotic procedure 
implies both 1/oj — > and 6 — > 0. It produces the multiplicity of available asymptotic 
paths on the plane (l/u,6); we parameterize these paths with a single small parameter 
e a = l/uj a , a = const > 0. This multiplicity gives us new opportunities for the study 
of a drift. At the end of the section we consider an important case when u does not 
contain any slow time-scale T: in this case we show the presence of an infinite number of 
solutions with a range of different T. 

In Sect. 3 we develop the inspection procedure for the governing equations. It leads to 
classification of different asymptotic solutions and to the key notions of critical, sub- 
critical, and super- critical asymptotic families of solutions. For a critical asymptotic 
family: (i) the small parameter is S\j2 (or a = 1/2); and (ii) the correlation (S3) between 
the oscillatory parts u and a of u and a enters the 0(1) approximation of the averaged 
governing equation. For a super-critical family a < 1/2 and again (S3) = 0(1), while 
some constrains on u must be imposed. For a sub-critical asymptotic family (ua) appears 
only in higher than 0(1) approximations. 

In Sect. 4 and Appendix A we obtain the general solution of an asymptotic problem for 
the critical asymptotic family (a = 1/2) with an arbitrary purely oscillating velocity u 
and arbitrary initial conditions for a. Each solution a = o+3 consists of explicit formulae 
for its oscillating part a and an advection-diffusion-type equation for the averaged part 
a. Our detailed analytical calculations include five steps of successive approximations. 
The obtained solutions satisfy the governing equations with small explicitly calculated 
residual right-hand-sides. It appears that evolution of the averaged field a includes a 
drift with the velocity V = Vq + eV \ + e 2 V '2 + 0(e 3 ) and pseudo-diffusion with a 
matrix of coefficients Ttjfc = £ 2 Xik + 0( e3 )i where e = £1/2 and Vo,Vi,V2,x itc have 
magnitudes O(l). Our term pseudo- diffusion means that the evolution of a is described 
by an advection-diffusion-type equation where the diffusion-type term appears not in the 
leading approximation. 

In Sect. 5 and Appendix B we present the asymptotic solutions for supercritical asymp- 
totic families with a = 1/3 and 1/4 and degenerated velocity field u. For a = 1/3 the re- 
quired degeneration is Vq = which leads to V = V\+eV2 + 0(e 2 ) and Kik = + O(e 2 ), 
where e = £1/3. For a = 1/4 the required degeneration is Vq = and Vi = which leads 
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to V = V2 + 0(e) and = + 0(e), where e = £1/4. Hence, the vanishing of the main 
term of a drift velocity Vq = (or the first two terms Vq = and V\ = 0) does not lead 
to the smallness of a dimensionlcss drift. Instead, a different expression V = 0(1) appears 
from a different asymptotic procedure. At the same time any supercritical oscillations do 
not produce pseudo-diffusion (within the precision of our calculations). 

In Sect. 6 we generalise the results of Sects. 4 and 5 considering more general u, which 
represents a power series of a small parameter e and contains both oscillating and non- 
oscillating parts. The expressions for drift velocities appear to be much more general, 
while pseudo-diffusion stays the same as in Sect. 4- 

Sect. 7 contains eleven examples of drift velocities and/or pseudo-diffusivity calculated 
for various flows. These examples are: (1) the velocity field u consisting of two arbitrary 
modulated fields of the same frequency; (2) u that represents the Fourier series of mod- 
ulated fields; (3) the Stokes drift, which allows us to make comparison (in Sect. 10) with 
classical results; (4) m is a spherical 'acoustic' wave; (5) Vi-drift (Sect. 5), which takes 
place for Vq = 0; (6) u is a plane travelling wave of a general shape or the superpo- 
sition of two such waves; (7) the polynomial velocity u, which has a nonzero average 
and represents polynomial of a small parameter; (8) u as the Bjorkncs configuration of 
two pulsating point sources; (9) a simple u that produces a drift Vq leading to chaotic 
dynamics of material particles; (10) the rigid-body-type oscillations of a fluid where a 
pseudo-diffusion term in the equations appears naturally; this example explains both 
mathematical and physical nature of pseudo-diffusion; (11) we consider a sub-critical 
asymptotic procedure and show that the averaged equations for a contain pseudo-drift 
and pseudo-diffusion; it emphasizes the significance of the critical and super-critical pro- 
cedures. 

The TTAM used in Sects. 3-7 combines five components: (i) a version of the two- 
timing method; (ii) the averaging over fast-time; (iii) the inspection procedure with the 
multiplicity of asymptotic procedures; (iv) the obtaining of general asymptotic solutions 
for these procedures; and (v) the estimating of RHS-rcsiduals (which has been only 
mentioned in Sects. 4-6 and is actually given in Appendices A,B). 

In Sect. 8 we study the links between the TT.AM-results of Sects. 2-7 and other theories: 
(i) the dynamical systems approach, see Aref (1984) Ottino (1989) , Samelson & Wiggins (2006) I; 



(ii) GLM-theory by|Andrews fc Mclntyre (1978)| (iii) the classical drift theory by |Stokes (1847)[ 



Longuet-Higgins (1953) Batchelor (1967)| (iv) the Krylov-Bogoliubo v averaging method, 

see |Bogoliubov fc Mitropolskii (1961)[|Krylov fc Bogoliubov (1947)||Sanders fc Verhulst (1985)[ 

(v) the homogenization method, see Bensoussan, Lions and Papanicolaou (1978) , Berdichevsky, Jikov, and Papanikolao 



(vi) the theory with nonzero molecular diffusion; and (vii) the MHD-kinematics in oscil- 



lating flows, see |Moffatt (1978)| The links between our TT^M-results of Sects. 2-7 and 
the theories (i), (ii) can be established only if we adapt and develop last two theories. 
Therefore, in Sect. 8.1 we calculate a drift velocity by the averaging of characteristics for 
transport equations and in Sect. 8. 2 we develop a two-timing version of GLM-kincmatics. 
It is important that the TTA M-thcory of Sects. 2-7 uses the Eulerian averaging while 
the theories (i) and (ii) operate with the Lagrangian averaging. These different averaging 
operations lead to different results for a drift and pseudo-diffusion. 



Finally, in Sect. 9 we discuss the assumptions used, the results obtained, and the prob- 
lems to be addressed. We also present short discussion/remarks at the end of each section 
and subsection. 
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2. Basic Equations and Operations 

2.1. Exact problem 

The equation for a scalar field a = a(x*,s*) in Cartesian coordinates x* = {x\,X2,x%) 
and time s* is 

JL + K.v.)) (*.,,.) = o (2 ,) 

where m* = u*(x* , s*) is a given velocity field, the asterisks mark dimensional variables. 
All functions used in this paper are considered to be sufficiently smooth. Equation (|2.ip 
describes the motion of a Lagrangian marker in either an incompressible or compressible 
fluid or the advection of a passive scalar admixture with concentration a(x*,s*) in an 
incompressible fluid. For a compressible fluid it also describes the advection of a passive 
scalar admixture, where a represents the ratio of concentration of admixture to density 
of a fluid. We study the motion of a Lagrangian marker (|2.1[) due to our primary interest 
in the motions of material particles. 

The velocity field u* has the functional form of a hat-junction 

u* = u*(x*,t*,T); t* = s* , t = u*s* (2.2) 

where t* and r are two mutually dependent time variables, which are introduced as two 
different time-scales k\s* and k2-s* with ki = 1 and ki = lu* (lu* is frequency). The r- 
depcndcncc in (|2.2I) is always 27r-periodic. The functional form (|2.2I) is aimed to describe 
modulated oscillations of high-frequency. Following established terminology we call t* 
slow time and r fast time. We do not require (|2.2[) to satisfy any equations of motion; 
such a general setting is often accepted, for example in the kinematic MHD-dynamo 
theory ( |Moffatt (1978) I. 

The functional structure of u* (|2.2[) underpins the idea that the solution of (|2.1[) also 
represents a hat-function: 

a = a(x*,t*,r) (2.3) 
Then after the use of the chain rule the equation (|2.1j) is 



_d_ 

ds* 



* d 




d 


+ (u 


a r 

OT 




d 






d 








+ 

i* 







3>*a=( w *_ + _ + (S*.V*))3 = (2.4) 

d d 



dr dt* 



(2.5) 



We also denote partial derivatives by subscripts, e.g. a T = da/dr. 

We obtain the dimensionless form of (|2.4p by introducing the characteristic scales of 
slow time T, length L, and velocity U. An example of a velocity field is 

u* = U(l + sin r)(l + C sin(t*/T)) g(x*/L, y'/i, (2.6) 

where g is an arbitrary vector-function of three scalar variables (one can also take T = 
L/U, but it narrows the class of velocity fields). Dimensionless variables and parameters 
are not asteriated: 

t = t* jT, x = x*/L, u = u*/U, lu = lu*T (2.7) 

According to Buckingham's 7r-theorem the problem (|2.4[) - (|2.7t possesses two independent 
dimensionless scaling parameters 

£! = l/u*T = l/u, 6 = U/(lu*L) (2.8) 



Admixture and Drift in Oscillating Fluid Flows 5 

where 5 can be interpreted as either of two physical parameters: (i) characteristic dimen- 
sionlcss displacement 5 = 5*/L (for displacement 6* = U/oj*); or (ii) Strouhal number 
5 = Q/uj* (for frequency fl = U/L). The parameter Si represents another Strouhal 
number and follows a more general notation: 

e a = l/u> a , with a = const > 0. (2.9) 

The dimcnsionlcss versions of (|2.5[) . (|2.4[) arc 

Da= [uj— + — +ujSu-V)a = (2.10) 
\ or at I 



d__ d_ 
ds 8t 



d_ 



Eqn. (j2.10[) can also be written in the form containing only small parameters: 

Da/uj = a T + e{a t + Su ■ Va = (2.12) 

One can see that we operate on the plane (e a , S) of two independent dimensionless scaling 
parameters. To use a rigorous asymptotic procedure one has to choose a one-parametric 
asymptotic path on this plane. In order to form such a path one might prescribe the 
dependence of each characteristic parameter in (|2.7[) on the chosen single parameter e a . 
The simplest choice is 

T = const, L = const, UT/L = 0(e a ) (2.13) 
We will exploit different options for paths, which go to the same limit 

(£ a ,S(e a )) -> as e a -> ( 2 -14) 
Remark: More details on the scaling procedure are given in Sect. 2. 3 and Sect. 7. 3. 

2.2. The classes o/H,B,T, and O(l) -functions 

Notations and definitions: 

Definition 1. The class H of hat- functions is defined as 

f{x,t,T) = f{x,t,T + 2it) (2.15) 

where t = s and r = los are two mutually dependent time variables (cj* is dimensionless 
frequency); the T-dependence is always 27r-periodic; the dependencies on x and t are not 
specified. 

Comments: (A) We have already accepted that u* e H (|2?2|) and a e H ((2~3| . (B) For 
any JeHa partial time-derivative can be expressed via the chain rule as 

— = ( — f(x,t,r) = f t + uf T , where t = s, r = us (2.16) 



9s \ (9i <9t 

(C) In any version of the two-timing method the variables t and r are considered to be 
mutually independent after the chain rule (|2.16[) has been applied. The justification of 
this auxiliary assumption is given a posteriori (see Appendix A). 

Definition 2. For an arbitrary / € H the Eulerian averaging operation is 



1 



TO+27T 



if) =2^ J /(*,*, r)dr, Vr (2.17) 
where during the integration t = const and (/) does not depend on tq. 
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Definition 3. The class T of tilde- functions is such that 

/ST: f(x,t,r) =f(x,t,T + 2ir), with (/) = 0, (2.18) 

Comments: (A) One can see that a T -function represents a special case of H-function 
(|2.15|) with zero average (|2.17j) . (B) In the text the tilde- functions are also called the 
purely oscillating functions. 

Definition 4- The class B of bar-functions is defined as 

/el: / r = 0, J(x,t) = (J(x,t)) (2.19) 

Comments: (A) Any M-function depends on x and t and does not depend on t; in 
particular, B-function can appear from H-function after its averaging. (B) With the use 
of the average (|2.17[) any H-function can be uniquely separated into its B- and T- parts: 

f(x,t,T)=J(x,t) + J(x,t,T) with (f(x,t,r)) =J(x,t) (2.20) 

Definition 5. The T -integration: For a given / we introduce a new function f T called the 
T- integral of /: 

f T = l f&'t'^-^l (£ J(x,t,o-)do-)d». (2.21) 

which represents the unique solution of the partial differential equation df T jdr = / with 
the condition (/) = (/ T ) = i|2T8]) . 

Comments: (A) The r-dcrivative of T-function always represents T-function. However 
the T-intcgration of T-function can produce an H-function. An example: / = f 1 shiT 
where f 1 is an arbitrary function of x,t: one can see that (/) = 0, however ( f£ fdr) = 
f 1 7^ 0, unless f 1 = 0. Formula (|2.21j) keeps the result of integration inside the T-class. 
(B) The T-integration is inverse to the r-diffcrentiation 

{f)r = (frV = /• (2.22) 

The proof is omitted. 

Definition 6. A dimensionless function / = f(x, t, r, e) (e is a small parameter) belongs 
to the class O(l) if / = O(l) and all partial x-, t-, and T-derivatives of / (required for 
further consideration) are also O(l). 

Comments: In all text below: (A) The highest spatial derivatives will be of the fourth 
order, while all t- and T-derivatives will be of the first order; (B) All large or small param- 
eters (in all this paper) are represented by various degrees of uj only; these parameters 
appear as explicit multipliers in all formulae containing tilde- and bar-functions; these 
functions always belong to 0(l)-class. 

Definition 7. The commutator of two vector fields a(x) and b(x) is 

[a,b} = (b-V)a- (a-V)b. (2.23) 

Comment: The commutator is antisymmetric and always satisfies Jacobi's identity for 
three vector fields a(x), b(x), and c(x): 

[a, b]=- [b, a], [a, [6, c]] + [c, [a, b}} + [b, [c, a]] = (2.24) 

Two more useful properties of commutators are: 

diva = 0, divb = => div[a,b]=0 (2.25) 
a-n = 0, bn = [a,b]-n = (2.26) 
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where n is a normal vector to a surface cW7. 
Properties of t -differentiation and T -integration: 
For r-dcrivativcs it is clear that 

fr =7r + fr = fr, (fr) = (fr) = (2.27) 

The product of two T-functions / and g forms a H-function: fg = F, say. Separating 
T-part F from F we write 

F = F-(F)=fg-(fg) = fg = {fg} (2.28) 

Since we do not use a two-level tilde notation for the T-parts of long expressions, we in- 
troduce braces instead. As the average operation (|2.17p is proportional to the integration 
over r, for products containing functions f,g, h and their derivatives we have 

( ? r ! = (( ? )r) _ = = ~ ( ^ t?> (2 - 29) 

{Ugh) + (fg r h) + (fgh T ) = (2.30) 
(f¥) = ((f T g T )r) - (f T g) = -(f T g) = -<f?) (2.31) 

which represent different versions of the integration by parts. Similarly 

([aX]) = -([a T ,b]) = -([5r,6]), ({aX}) = -([a T ,b]) = -({a T ,b}) (2.32) 

The equation to be solved: 
Let us accept that a dimcnsionless velocity in (|2.10p . (|2.12l) is 

u(x,t,r) einO(l), and 5(ui) = Squj^ 1 , /3 = const < 1 (2.33) 

where u is a given function; indefinite constant /3 allows us to consider different asymp- 
totic paths; one can take 5o = 1 since any other <5 can be incorporated into u. Eqn. 
(|2.10[) can be rewritten as 

Da = uja T + a t +lj^(u ■ V)o = 0. (2.34) 

In Sects. 3-7 we will study and solve this equation. 

Remark: The r-periodicity of the given velocity (|2.2[) , (|2.33[) represents a restriction 
which is artificially imposed to simplify further calculations. At the same time we will 
show that a r-periodic asymptotic solution (|2.3[) can be build for any initial data; hence 
we are not forced to go outside the H-class (|2.15l) . 

2.3. Choice of scale T for t-independent velocity 
A i-indcpcndent velocity 

u = u(x, t) (2.35) 

is important theoretically. This case represents a degeneration, where the presence of 
slow time t in eqns. p,10|) . (|2.34p is provided only by the additional agreement a t £ O(l) 
leading to two mutually dependent time-scales (t, t) = (s,us). It is clear that for (|2.35|) 
this agreement represents only one available option and an infinite number of different 
slow-time-scales can be introduced. To make it possible one can write 

t = uj 1 ~ x (uj x s) = LO\t\ = t\ where u>\ = w 1_A , t\ = w A s, A = const < 1 — Ao(2.36) 

with a small constant Ao > which provides uj\ — > oo when u — > oo. These manipulations 
allow to introduce time-scales (t\,T\) = (u x s,ijjs) and a new version of eqn. (|2. 341) 

— +^ux-V)a = 0, — =uj x —+lj— (2.37) 

OS J OS Ot\ OT\ 
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with u\ = u(x,t\). Further transformations of (|2.37[) yield 

(J^+Ux^)a + u^(u x -V)a = 0, (2.38) 

where (3\ = ((3 — A)/(l — A) and a = a(x 7 t\ 7 T\). Finally, we replace the old scaling 
agreement at € ©(1) with a new one at x £ 0(1) (a Tx <E O(l) remains valid since t\ = t). 
Now one can see that the problem (|2.38[) is mathematically identical to (|2.34[) . which 
means that for any solution in variables (t, r) we have an additional solution in variables 
(t\,T\) (and vice versa); the varying of A gives us an infinite number of such additional 
solutions. Hence, an infinite number of new solutions to (|2.38p can be obtained by rescal- 
ing of the known solutions to (|2.34p . which will be considered in Sects. 3-7. The existence 
of such additional solutions means that motions with any slow-time-scale Lo x t (|2.36l) arc 
possible. Indeed, it sounds physically reasonable: if a slow-time-scale is not enforced ex- 
ternally (does not appear in the given velocity) then any time-scale that is longer than 
imposed oscillations can be considered as a slow one. For brevity, in the text below we 
keep at <£ O(l) and two time-scales (i, r) = (s,uis). 

Remark: Similar to (|2.35|) . (|2.38j) . the multiplicity of slow-time-scales takes place, for 
instance, in dynamics of a classical pendulum with a harmonically vibrating pivot, when 
an external gravity field is absent. 



3. The Inspection Procedure 

For simplicity we consider now cqns. (|2.33l) . (|2.34[) with purely oscillatory velocity u = 
u G TnO(l). In our inspection procedure we use a test solution to (|2.34l) 

a(x,t, t) = a(x,t) -\ — — b(x,t,r); a = const > (3.1) 

where nelfl O(l), b e T n O(l) and the amplitude of an oscillating part is given by 
the small parameter e a = l/tu a (|2.9j) . The substitution of (|3.1|) into (|2.34|) yields 

a t + — (ujb T + b t ) +UJ 13 ( (5-V)o+ — (u-V)6) =0 (3.2) 

LO a \ I \ UJ a J 

The B-part of (J3HJ) is 

a t +oj l3 - a ((u-W)b) =0 (3.3) 

and the subtraction of (|3.3[) from ()3.2|) produces the T-part of p.2j) 

uj^br + uj- a 1 t + ujP(u- V)q + uj p - a {{u ■ V)6} = (3.4) 

with the notation {•} (|2.28[) used. In each underlined group in (|3.4p the first term domi- 
nates over the second one (due to a > (|3.1[1 ). Therefore, to make this equation mean- 
ingful we need to accept that the dominating terms are of the same order: 

l-a = p or a + (3 = 1. (3.5) 

Then the imposed restrictions a > (|3.1j) and j3 < 1 (|2.33|) are mutually compatible and 
(|3.3j) . (|3.4|) can be rewritten as 

at + w 7 ((m • V)b) = where j = 2/3 -1 = 1- 2a (3.6) 
uj (X + (u ■ \ +b t + ^{{u ■ V)6} = (3.7) 
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The second term in eqn. (|3.6|) represents the changes of a due to the averaged nonlinear 
effects of oscillations. Taking the different values of 7 we introduce three qualitatively 



different asymptotic families of solutions: 

7 > (1/2 < j3 < 1, < a < 1/2) —super- critical oscillations (3-8) 

7 = (j3 = 1/2, a = 1/2) —critical oscillations (3-9) 

7 < (—00 < P < 1/2, 1/2 < a < 00) —sub-critical oscillations (3.10) 



One can notice that super-critical, critical, and sub-critical oscillations have been intro- 
duced as the asymptotic families of solutions, not as particular solutions (recall that any 
particular solution can be treated as a member of several different asymptotic families). 

For the family of critical oscillations the averaged oscillatory term in p.6[) has the 
same order O(l) as a mean solution. Then the leading terms (for large w) in (|3.6[) . (|3.7p 
are: 

3t + <(5-V)6)=0; (3.11) 
& r + (w-V)a = 0. (3.12) 

This system gives a single equation for a 

Ot-((«-V)(C-V)>o = 0, £ = u T (3.13) 

which is obtained by T-integration (|2.21[) of (|3.12j) and the following substitution of 
b(x,t,r) = -(£ ■ V) a into (j3"TT]) . 

The critical family of oscillations (|3.9[) separates the larger families of sub-critical and 
super-critical oscillations from each other. For super-critical oscillations 7 > and the 
correlation in (|3.6p can not be balanced by a t = 0(1). Then the condition of solvability 
requires this correlation to vanish, hence the leading order of (|3.1ip . (|3.13[) is replaced 
with {(u • V) 6) = 0. It means that the imposed velocity must degenerate (satisfy certain 
additional restrictions) . 

Remarks: 

1. The inspection procedure is aimed to find and classify all possible asymptotic solu- 
tions (of the functional form (|3.1[) ) of a given PDE. This procedure also allows to choose 
a natural small parameter (|2.9p with different values of a for different classes of solutions. 

2. The critical and super-critical oscillations are most interesting both theoretically 
and practically. We justify this statement with three arguments: (A) Sects. 4-7 below are 
completely devoted to these two oscillatory families. The obtained results are mathemat- 
ically satisfactory and look physically convincing. (B) These families of oscillations are 
the least studied, probably due to the reason that for them U — > 00 as uj — > 00 in f)2 . 13[) ; 
however it is well-known that such an increase of U does not diminish the validity of 
asymptotic procedures. (C) In Sect. 7. 11 we will consider the sub-critical family (|3.10p 
with f3 — in ()2.33[> . (|2.34|) . We will show that for this family a drift is replaced with a 
pseudo-drift and the general solution is diverging. 



4. Family of Critical Oscillations. Purely Oscillatory Velocity. 

4.1. Critical asymptotic procedure with a = 1/2 
Following (|3.9p we accept that f)2 .33[) belongs to a critical asymptotic family: 

u(x,t,r) = u(x,t,r), 5 = u>~ 1/2 , (4.1) 
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Then (|2Jll) , 63) yield: 

2>a = a;a r + at + \fu{u ■ V) a = (4.2) 

The small parameter £1/2 = l/v^ (|2.9|) allows to rewrite it as 

D 2 a = a T + e{u- V)a + e 2 a t = 0, £> 2 = £ 2 D = D/w (4.3) 

where the subscript in Ei/ 2 has been dropped. We are looking for the solution of (|4.3[) in 
the form of regular series 

00 

a(x,t,r) = ^2e k a k (x 1 t,r), 2 t elnO(l), ft = 0,1,2,... (4.4) 

fc=0 

The substitution of (|4.4p into (|4.3I) produces the equations of successive approximations 

a 0T = (4.5) 

ai T = -(u. V)2o (4.6) 

a„ T = -(«• V)a„_i - d t a n -2, d t = d/dt, n= 2,3,4,... (4.7) 
S fe = a fe (x,t) +a fc (x,t,r), a fc G BnO(l), S fc G TnO(l), ft = 0,1,2,... 

where (|4.6p represents a systematically derived version of ([3~T2|) . Let oW be a truncated 
solution (|4.4[) of order n: 

n 

a w (a;,t,r) = ^ e fc a fc (a;, t, r) (4.8) 

fe=0 

Its substitution into (|4.3|) produces the RHS-rcsidual Rcs[n]: 

D 2 al"] = Res[n] (4.9) 

4.2. General solution for the first five approximations 

The detailed solving of (|4.5[) - (|4.7j) for n = 0, 1, 2, 3, 4 is given in Appendix A. The trun- 
cated general solution is 

= a + eoi + £ 2 a 2 + £ 3 a 3 + £ 4 a 4 (4-10) 
with the tilde-parts given by explicit recurrent expressions 

5 = 0, (4.11) 

51 = -(€-V)ab, (4.12) 

5 2 = -(rV)a 1 -{(i t -V)S 1 } r , (4.13) 
a 3 = -(i-V)a 2 -{(«- V)S 2 } T -Sl t , (4.14) 
04 = - (€ ■ V) S3 - {(« • V) a 3 } T - a T 2t , (4.15) 

£ = u T = [ u(x,t,a)da — — [ ( [ u(x,t,a)da\ dfi (4.16) 
Jo 2tt J \J J 
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and the bar-parts au satisfy the equations 

(d t + V o -V)a o = (4.17) 
(d t + V -V)ai + (Vi ■ V)a = (4.18) 

(8 t + Vo ■ V) a 2 + (Vt ■ V)aj + (V 2 • V)a = ^ (**|^) ' ( 4 ' 19 ) 

V = I[u,e], V 1 = \([[u,lU\), (4.20) 

F 2 = \([[V Q ,lU]) + i([Vo, Vol) + 5<K»€t]> + i«divu' + S'divO, (4.21) 
5' = | t -[V 0j ^], (4.22) 
2x ifc = + u'fcfi) = £ Fo , (4.23) 

% 7« = (ft + V„ • V) - f^7 m - ^T-lkm (4-24) 

ox m ox m 

where the operator £y is such that Suy fik = represents the condition for tensorial 
field fih( x i t) to be 'frozen' into velocity field Vq(x, t). Three equations (|4.17|) - (|4.19|) can 
be written as a single advection-pseudo-diffusion equation (with the error 0(e 3 )) 

^ +v ^'l,(^S) (4 ' 25) 

V = V l2] =V + eV x + e 2 V 2l (4.26) 

*** = Xffc = z 2 X lk (4.27) 
a = a [2] = a + eai + e 2 a 2 (4.28) 

where all coefficients are given in (|4.20[) - ([4.24[) . Eqn. (|4.25p shows that the averaged 
motion of a represents a drift with velocity V = V + 0(e 3 ) and pseudo-diffusion with 
the matrix coefficient Kn~ = e 2 \ik + 0(e 3 ). 

Definition: The term pseudo- diffusion (PD) means that: (i) the evolution of a is de- 
scribed by an advection-diffusion-type equation (|4.25j) : (ii) the diffusion- type term ap- 
pears in (|4.25j) not in the leading approximation (see (|4.19[) , where it represents a known 
source-type-term in the second approximation) ; and (ii) the equation (|4.25|) is valid only 
for regular asymptotic expansions (|4.26p - (|4.28p . 

In a' 4 ' (|4.10p the expressions for a , a\, a 2 , Si, a 2 , a 3 , and 04 are given by (|4.11|) - (|4.24|) . 
while 0,3 and 04 are to be found from higher approximations. The substitution of into 
(|4.3p produces the RHS-rcsidual of order e 5 

D 2 a [4] = Rcs[4] = 0(e 5 ) = 0(l/w 5/2 ) (4.29) 

One can notice that the residual of order 0(l/w 5 / 2 ) for D 2 leads to the residual 0(1/cj 3 / 2 ) 
for the original operator £) (|4.2p . 
Remarks: 

1. The statement that ()4. 1 ![) - f|4.24|) describe the solution of the problem with arbitrary 
initial data is based on the fact that initial data for (|4.17p - (|4.19p can be chosen arbitrarily. 

2. We have explicitly solved the equations of the first five approximations. Such persis- 
tence does not occur often in the use of asymptotic methods. However, the well-known 
examples are the surface gravity wave, see Stokes (1847)| Stoker (1957) Debnath (1994) 



and the low-i?e-solutions for a sphere, see Chester, Breach, and Proudmen (1969) Our 
motivation for doing the cumbersome calculations of Appendix A is based on two our 
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results: (i) the first and the second approximations of the critical asymptotic solutions 
appear as the main-order terms in super-critical solutions (see Sect. 5); and (ii) PD (which 
qualitatively complements a drift) appears only in the second approximation. 

3. The term PD is different from super- diffusion, which means the spreading of an ad- 
mixture faster than \ft (see Avallaneda &: Majda (1992)| Volpert, Nec, and Nepomnyashchy (2010) |. 



5. Super-critical Oscillations 

5.1. Super- critical asymptotic procedure with a. = 1/3 

For critical oscillations (j4~Tj) the drift velocity (|426|) is V ~ V = 0(1) P~2"0]l . From 
(|4.26[) one can expect that for the degenerated u such that 

P f o = i([S,€]) = (5.1) 

the drift is V ~ eVi = 0(\/^/uj). However, one can find that the drift velocity is 
still 0(1) in a different asymptotic family. To show it, we accept (|5.1|) and consider a 
super-critical family of purely oscillating velocities 

u(x, t, r) = u(x, t, r), 5 = uj~ 1/3 (5.2) 

where u(x,t,r) 6 Tfl O(l) is a given function. Then eqns. (|2.34l) . (|5.2[) yield: 

Da = uia T +a t + oj 2/3 (u ■ V)a = (5.3) 

A small parameter ei/3 = cj -1 / 3 allows to rewrite ()5.3[) as 

T) 3 a = a T + e(u ■ V)a + e 3 a t = 0, D 3 = e 3 D = D /lj (5.4) 

where the subscript in £1/3 is dropped. We are looking for a solution of (|5.4j) in the form 
of regular series (|4.4[) with the redefined e. The substitution of (|4.4[) into (|5.4|) produces 
the equations for successive approximations (cf. with (|4.5[) - (|4.7[0 

aor = (5-5) 

a\r = — (u ■ V) ao (5-6) 

a 2r = -(2 • V) ai (5.7) 

a n r = -(« • V) a n _i - <9 t a„_ 3 , n = 3,4,5,... (5.8) 

These equations can be solved in the most general form (see Appendix B). Here we 
present the truncated solution a' 4 ' (|4.10p . For the tilde-part of each approximation one 
can obtain explicit recurrent formulae 

50 = 0, (5.9) 

51 = -(|-V)a , (5.10) 
a 2 = -(€ -V)oi -{(«■ V)S!} T (5.11) 

53 = -(£■ V)a 2 -{(«■ V)S 2 } T (5.12) 

5 4 = -(€ ■ V)o 8 - {(« • V)S 3 } r - a[ t ,Va 3 (5.13) 
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and the transport equation (which is valid with the error 0(e 3 )) for averaged motion 

(d t + V ■V)a = 0, (5.14) 

V = V [2] =V 1 +eV 2 , (5.15) 
a = aP^ =a a + eai, Vc^ 03,04 (5.16) 

V! = -±{[K,t]), V 2 = - 1 -{Tl + K n ), (5.17) 

X= [£,«], k= [K T ,K], K" = [[K,$U] (5.18) 

Eqn. (|5.14[> means that the averaged motion of a passive admixture (up to the first two 

[2] 

approximations) represents a pure drift with velocity V (|5.15|) . 

Hence, in (|4.10p functions aoi oa, 02, S3, and 04 are given by (|5.14[) - (|5.17D 
and l|5.10p - (|5.13p . while 0,2, 03 and 04 are to be found from further approximations. The 
substitution of into (|5.4p produces the RHS-residual of order e 5 

33 3 aW = R CS [4] = (9(e 5 ) = 0(l/a; 5/3 ) (5.19) 

The residual of order <3(l/u; 5 ' 3 ) for the operator D3 means that the residual for original 
operators ([23]), (123), (12331) . and flJTS) is 0(l/w 2 / 3 ). 

5.2. Super- critical asymptotic procedure with a = 1/4 

Next we consider the oscillating velocity fields with first two terms of the drift velocity 
1)4.261) vanishing 

V = = and y 1 = I([[u,i],|])=0 (5.20) 

Here one can find that the drift of magnitude 0(1) appears within the asymptotic family 
with a = 1/4. Let us accept (|5.20p and introduce a family of super-critical oscillations 

u(x, t, t) = u(x, t, r), 5 = oj- 1/4 (5.21) 

where u{x, t, r) e T n O(l) is a given function. Eqns. (|234]) . ([5T2"T]l give: 

Da=wa T +fl t + w 3/4 («'V)a = (5.22) 

The small parameter E\u = uj^ 1 ^ 4 allows us to rewrite it as 

£> 4 a = a T +e(M- \7)a + e 4 a t = 0, £4 = £ 4 £> = £>/w (5.23) 

where the subscript in ex/4 is dropped. We are looking for a solution of f|5 -23[) in the form 
of regular series f|4.4[) with the redefined e. The substitution of f|4.4[) into (]5.23l) produces 
the equations for successive approximations (cf. with (|4.5p - (|4.7p and (|5.5p - (j5.8p ) 

a 0T = (5.24) 

a lT = — (tt-V)a (5.25) 

a 2r = — (u-V)ai (5.26) 

a 3T = -(M-V)a 2 (5.27) 

a„ T = —(u ■ V) a n _i - <9 t a„_4, n = 4,5,... (5.28) 

Again, these equations can be solved in the most general form (see Appendix B). Here 
we present the truncated solution a" (|4.10p . For the tilde-part of each approximation 
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one can obtain explicit recurrent formulae 

S = 0, (5.29) 

ai = -(€-V)5o, (5.30) 

5 2 = -(£■ V)oi -{(«■ V)Si} r (5.31) 

5 3 = -(€ ' V)a 2 -{(«■ V)S 2 } r (5.32) 

5 4 = -(£-V)a 3 -{(S-V)S 3 } r , Va 3 (5.33) 

and the transport equation for the averaged motion 

(d t + V -V)a = 0, (5.34) 
V = V m = V 2 , a= a [2] = a ; 



with the same notations as in (|5.17[) . (|5.18|) . Equation (|5.34|) means that averaged motion 
represents a pure drift with velocity V 2 . 

Hence, aM can be written as (|4.10p where So, Si, So, S 2 , S 3 , and S4 are given by 
(|5.29p - (|5.34p . while 01, a 2 , a 3 and 04 are to be found from further approximations. The 
substitution of a' 4 ' into (|5.23[) produces the RHS-residual of order e 5 

D 4 SW = Rcs[4] = 0(e 5 ) = 0(1/lu 5 / 4 ) (5.35) 

The residual of order 0(l/w 5 / 4 ) for the operator Q± means that the residual for original 
operators ([23]) . CO}, (|2~33|> , and ([5T2"2"|) is (^(l/w 1 / 4 ). 

Remark: One can see that the results for both the super-critical families are simpler 
than for the critical one (Sect. 4) and they are build up from the same elements: if the 
conditions of degenerations (|5.1|) or (|5.20[) are valid, then either V\ or V 2 (|4.26|) becomes 
the leading 0(l)-term in (|5.16|) or (|5.35j) . At the same time pseudo-diffusion does not 
appear (within the considered precision). 



6. Polynomial Velocity for Critical and Super-Critical Oscillations. 

Expressions (|4.1[) , (|5.2I) . (|5.2ip for velocity have a special form. One can expect that 
the adding of higher-order terms to these expressions may potentially change both the 
drift and pseudo-diffusion. 

6.1. Critical asymptotic procedure with a = 1/2 

Let us consider a critical asymptotic family with a velocity in (|2.2p . (|2.10p . (|2.34p given as 
(cf. with 03])) 

u = v H — j=w -\ — r. S = uj~ 1 / 2 1 (6-1) 

where v,w,f £ H D O(l) are given functions. This formula might be seen as the first 
three terms of an infinite series, truncated by the omitting of terms 0(l/w 3 / 2 ). Eqns. 
(EM)), (EH) yield: 

Da = luo, t + y/uj(v ■ V) a + (d t + w ■ V)a + — = (r ■ V) a = (6.2) 

The small parameter Syi 2 = l/\/w allows to rewrite it as 

D 2 a = a T +e(v-V)a + e 2 (d t + w ■ V)S + e 3 {r ■ V) a = (6.3) 
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where D2 = s 2 T) = T>/u>, and the subscript in £1/2 is dropped. We are looking for the 
solution of (|6.3|) in the form of regular series (|4.4j) . The substitution of (|4.4j) into (|6.3[) 
produces the equations for successive approximations 

a r = (6.4) 

air = — (v ■ V) ao (6.5) 

a 2T = — (v ■ V)ai — (d t + w ■ V)a (6.6) 

a.3r = — (v ■ V)a 2 — (dt + w ■ V)ai - (f • V)ao (6.7) 

The solving of (|6.4p - (|6.7[) follows the same steps as in Appendix A. An important new 
element is: the bar-part of eqn. (|6.5[) is v ■ Vao = 0. It requires v = as a solvability 
condition. Therefore the leading term in the prescribed velocity (|6.1|) must always be 
purely oscillatory. 

Here we present only the truncated solution a' 3 ' (|4.10p . For the tilde-parts with 
k = 0, 1, 2, 3 the equations (|6.4p - (|6.7[) yield the recurrent expressions 

50 = 0, (6.8) 

51 = -(i-V)a , (6.9) 
a 2 = -{(v ■ W)ai} T - (£■ V)ai - (t? ■ V) a (6.10) 
S3 = -{(« ■ V)S 2 } r - {(«> • W)ai} T - (6.11) 
-(C " V) a 2 - (n ■ V) 01 - (C ■ V) a - (ft + ^ ■ VJoJ 

I = w r , n = w T , C = r T , (6.12) 

and for the bar-parts Ofc we obtain the transport equation 

(d t + (w + Vo)-V)a o = (6.13) 
(d t + (137 + F ) ■ V) Si + ((r + Vi + V12) ■ V)a = (6.14) 

V = l([5,t]), Vi = l([[v,£U}), Vi 2 ={[v,?j\), (6.15) 
which can be rewritten (with the error 0(e 2 )) as a single advection equation 

(d t + V -V)a = (6.16) 
V = (w + V ) + e(r + Vi + Vi 2 ), a = a + eai (6.17) 

Notice that the solution given by (|6.8j) - (|6.17j) satisfies eqn. (|6.3p with Res[3] = 0(lo~ 2 ). 

Remark: In (|6 .4[) - (|6 . 1 7[) we present only four steps (k = 0, 1, 2, 3) of successive approx- 
imations since that is enough to give us an instructive correction V12 to the previously 
calculated drift ()4.18j) . (|4.20p . The fifth step (fc = 4) has been also carried out but for 
brevity it is not shown here: it produces a combination of a drift term V 2 (which is 
more cumbersome than (|4.2ip ) and pseudo-diffusion with the same coefficients (|4.23l) . 
The related calculations arc similar to those of Appendix A. 

6.2. Super- critical asymptotic procedure with a = 1/3 
Next, we consider a super-critical asymptotic family with a velocity similar to (|6.ip 

u = v + uj~ 1/3 w + w~ 2/3 r, 6 = uj~ 1/3 (6.18) 

where v,w,r G HnO(l) are given functions. Expression (|6.18|) might be seen as the first 
three terms in an infinite series, truncated by omitting the terms 0(l/w). Eqns. (|2.34p . 
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f6~T8| yield: 

£>2 = ua T + oj 2/3 (v ■ W)a + uj 1/3 (w ■ V)a + (d t + r ■ V)a = (6.19) 

The small parameter = w^ 1 / 3 (|2.9p allows to rewrite (|6 . as 

£ 3 a = a T + e(v ■ V) a + e 2 (u> • V) a + £ 3 (<9 t + f • V)a = (6.20) 

where D3 = e 3 D = D/u and the subscript in ei/3 is dropped. We are looking for the 
solution of (|6.20[) in the form of regular series (|4~4|) . The substitution of (|4~4)l into f|6.20[) 
produces the equations for successive approximations 

a 0T = (6.21) 

a lT = -(v-V)a (6.22) 

a 2r = — (« ■ V)ai — (u) ■ V)a (6.23) 

a 3r = -(v- V)a 2 - (w ■ V)ai - (9 t + f • V)a (6.24) 

The solving of (|6.21l) - (|6.24l) follows the same steps as in Appendix B. There are two 
solvability conditions 

v = 0, V = -w (6.25) 

which are required for the existence of the bar-part solutions in (|6.22|) and ()6.23|) . One 
might see the consideration of Sect. 5 as being physically incomplete. Indeed, one can 
expect that the results for small Vq 7^ are close to that for Vq = 0. This expectation 
has been met and fulfilled by a more general degeneration condition Vq = —w (|6.25j) 
which leads to a more general (than (|5 . 1 5[) ) expression for a drift. 

For tilde-parts one can obtain the explicit recurrent expressions which coincide with 
(|6.9[) - (|6.1ip . while the formula (|6.12[) for a 3 should be replaced by 

a 3 - -{(v ■ V) a 2 } T - {(w ■ V)Si} T - (6.26) 

-(£ • V)a 2 - (?7- V) ai - (?• V)a - (to • V)oJ 

For the bar-parts one can derive the equation 

(dt + V ■ V)a = 0, V = r + Vi + V 12 (6.27) 

which is valid with the error 0(e 2 ) and contains the same Vi and V12 as in f)6. 15[) . 
Notice that this solution satisfies the equation (|6.20[) with Res[3] = 0(w -4 / 3 ). 
Remarks: 

1. One might additionally consider a polynomial version of velocity for w^-family 
(|5.21| . The related results could be predicted by the analogy between the material of 
Sect. 6. 2 and Sect. 5. 

2. For the critical family (|6.1[) slow motion represents the drift (|6.17p of magnitude 
0(1) combined with the same pseudo-diffusion of order 0(l/u>) as in (|4.28p . hence the 
results are very similar to those of Sect. 4- 

3. All the results of Sects. 2-6 can be generalized to an arbitrary flow domain f2 with 
fixed boundary dfl. Then (|2.26p yields that 

V ■ n = Vi ■ n = V 2 ■ n = u ■ n = on dQ (6.28) 

provided u satisfies the no-leak condition. One can also see that if a fluid is incompressible 
then (|2.25p yields that all drift velocities in (|4.17|) - (|4.24|) are also incompressible: 

divw = => div V = div V x = divF 2 = (6.29) 

The explicit expressions (|4.17|l - (|4.22|) show that due to incomprcssibility only the last 
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term in V2 (|4.21l) vanishes. Similar (for (|6.13[) - (|6.15[0 the incomprcssibility yields 

divv = divw = divr = =>■ div Vq = div V\ = divVi2 = (6.30) 

4. For the super-critical family (|6 . 18[) with the degeneration (|6.25[) slow motion is a pure 
drift with the velocity (|6.27[) which is 'upgraded' to O(l) from 0(e)-term in (|6.17|) . The 
qualitatively new addition to a drift (in comparison with a critical family) is a cross-term 
V\2 = ([v, rj]). It represents an averaged commutator between two mutually independent 
functions v and rj; hence V12 expands the class of available analytical expressions for a 
drift. For example, for an incompressible fluid ([v, rj]) = V x (v x rj), so one can see that 
(v x rj) represents an arbitrary vector-potential for V12. 



7. Examples. 

7.1. Superposition of two modulated oscillatory fields of the same frequency 
The velocity field u gT]) and | = u T ([2~2T]) are 

u(x,t,T) =p(x,t) sin t + q(x,t) cost (7-1) 
£(x,t,r) = —p{x, t) cos r + q(x, t) sin r (7.2) 

with arbitrary B-functions p and q. The straightforward calculations yield 

[«,£] = M (7-3) 
so the commutator surprisingly is not oscillating. The drift velocities (|4.20[) . (|4.21[) are 

Vo = l([u,^}) = \[p,q}, Vi = \{[[n,M]) = 0, (7-4) 
V 2 = 1 ([P,p] + [Q,q]) - \ ([p t ,p] + [q t ,q]) + (7.5) 

+ ^ (p div T' + qdiv Q' + P'div p + Q'div q^j , 

P=[V Ql p], Q = [V ,q], P'=p t P, Q' = q t Q, 

(Mk) = -^{PiPk+lilk) (7-6) 

A pseudo-diffusion matrix , which follows after the substitution of (|7.6p into (|4.23l) . 
Remarks: 

1. The degeneration Vq = in an incompressible fluid (|6.29[) , (|6.30p corresponds to 
p xq = with an arbitrary potential Tp(x, t). For any such fields all results of Sect. 5.1 
are valid. 

2. The velocity u (|7.ip is general enough to produce any given function Vo(x,t). To 
obtain p(x,t) and q(x,t) one has to solve a bilinear first-order PDE 

[p, q] = (q ■ V)p - (p ■ V)q = 2V (x, t) (7.7) 

which represents an undcrdctcrmincd bi-lincar PDE-problcm for two unknown functions. 
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7.2. Fourier series of modulated fields 
The previous example can be generalized by the consideration of a velocity 



u(x, t, t) = Pk( x i t) sm ^ T + Qk( x i t) cos kr (7-8) 
fe=i 

CO 

fe=i 

The calculations of Vo lead to an infinite sum of commutators 

1 00 1 

^o = 2<M>=E2fc[M fc ] (7-10) 
fe=i 

In particular, (|7.10j) gives an infinite number of fields u with Vo = 0. For example, if 
p k = 0, Vfc then 

m±n±/— 

where the sum is taken over all sets of three positive integers m, n, I such that m±n±Z = 0. 
It is clear that for many combinations Vq = and V\ ^ (see Sect. 5.1). 

7.3. The Stokes drift 

This most celebrated example of a drift we consider in some details. We start with the 
scaling ()2.7|) . (|2.8[) and the limit (|2.14j) . The dimensional solution for a plane potential 
travelling gravity wave (see |Stokes (1847)|[Lamb (1932)] |Debnath (1994)[ ) is 

S . =M , a = e x p( W (^; : ;)) (,12, 

where ui*, k* , /i*, and g* are dimensional frequency, wavenumbcr, spatial amplitude, and 
gravity; (x*,y*) are Cartesian coordinates, and r = to*t* . One can immediately notice 
that: (i) the characteristic length is L — 1/k*; and (ii) the scale U is apparent. Hence the 
dimensionlcss scaling parameters (|2.8[) . (|2.12[) and the asymptotic limot (|2.14[) appear as 
k = 1 and 

E X = 1/Tu* = l/w 0, S = U/cv*L = g*h*/L 2 uj* 2 = gh/uj 2 -> as w oo (7.13) 

where one can choose the scale T = T(lu*) in an arbitrary way, its only mission is to 
provide E\ — > as ui — > oo (see Sect. 2. 3). The simplest way to provide S = (|4.1[) is 

to choose an asymptotic family as 

T = const, L = const, gh = 0{uj 3/2 ) (7.14) 
The dimensionless velocity field (|7. 1 2[) and £ are 

u = Ae ky ( °° s( - kx ~ T M £ = ^ ( ~ Sin(fc ' T ~ T M (7 15) 

Y sui(kx — r) y ' I cos(fcx — r) y 

where in the chosen system of units A = 1 and fc = 1; however, we keep both A and A: in 
the formulae for tracking its physical meaning. The fields p(x, y), q(x,y) (|7.1|) arc 

— a kii ( sin kx \ _ . ki, f cos fa \ m , „\ 

P = -4e fcy . , q = Ae fc!/ . , (7.16) 

i — ™c kx V sin kx ' 



COS KX 
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= fcA 2 e 2fcy ylj, Vi = (7.17) 

which gives Vq proportional to the classical Stokes drift and a zero value for the first 
correction to it; the explicit formula for V2 is not given here for brevity. At the same 
time 

(tek)=Z{x,y,t)6 ik , with E=^A 2 e 2ky (7.18) 

which can bee seen as 'locally isotropic' oscillations. The matrix of pseudo-diffusion (|4.23p 
is 

dV i , dV k 



2 Xik = { (dt + V ■ V) S ik - 2x e ik } 5, 2e ik = 1- 
Further calculations show that 



__ / 1 



with x= -k 2 A i e 3ky 
10/ 4 



One can see that the eigenvalues X\ = — X an d X2 = X correspond to ordinary diffusion 
in one direction and anti-diffusion in the perpendicular direction, so here we have a 
mixed case of pseudo-diffusion. The averaged equation (|4.25[) (with an error 0(e 3 )) can 
be written as 

a t + (Vo + e 2 V 2 )a x = e 2 (x y a x + xa xy ) (7.19) 
a = clq + ea\ + e 2 a,2 

where Vq and V2 are the x-components of corresponding velocities (their y-components 
vanish). This equation has an exact solution a = a(y) (with an arbitrary function a(y)), 
which is not blurred by pseudo-diffusion. 
Remarks: 

1. In order to avoid confusion one should recall that the choice of asymptotic family 
(|7.14[) has nothing to do with physical relations between parameters. Any asymptotic 
family represents only a formal one-dimensional parametrization in the functional space 
of all possible functions u* (|7.12|) . The only aim of such a parametrization is to build 
a valid asymptotic procedure. Therefore any particular wave (|7.12[) with given values 
k* = /cq and u* = Wg can be considered as a point on an asymptotic curve with k = 1 
and with variable u irrespectively to the presence of any dispersion relation between fcg 
and coq. 

2. A parameter gh — > 00 as u> — > 00 (|7.14p which is a common property for many 
asymptotic procedures. 

3. More general (than ()7.14[) ) asymptotic families are considered in Sect. 2. 3 and in 

4. Both fields (|7.15[) arc unbounded as y — > 00, but it is not essential for our purposes. 

7.4. Spherical 'acoustic' wave 
A velocity potential for an outgoing spherical wave is 

0=-sin(fcr-T) (7.20) 
r 
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where A, k, and r are an amplitude, a wavenumber, and a radius in a spherical coordinate 
system. The velocity is purely radial has a form ()7,ljl 



u = psinr + gcosr, (7-21) 
p = A ^— cos kr H — sin fcr^ , 11= A ( — - sin fcr H — cos kr^j (7.22) 

where u, p, and q are radial components of corresponding vector-fields. The fields £ and 
[u, £] arc also purely radial; the radial component for the commutator is 

£u r - u£ r = A 2 k 3 /r 2 (7.23) 

where £ is radial component of £ and subscript r stands for radial derivative. The drift 
(f773|) . (f7T4|) is purely radial with 

^o = 4^, F 1= 0, V 2 = ^(3k>-^ (7.24) 

It is interesting that Vq formally coincides with the velocity caused by a point source 
in an incompressible fluid and for small r the convergence of V is endangered, since V2 
dominates over Vq. Further calculations yield 

<£ 2 > = ^( fc2 + 1 A 2 ), X = A 4 k 5 /4r 2 >0 (7.25) 

where x stands for the only nonzero rr-component of Xik ■ One can see that in this case 
pseudo-diffusion appears as ordinary diffusion. 

7.5. TheVi-drift. 

If Vo = then the drift of order O(l) is given by ([5715]) . Let the velocity field d572j) be 
a superposition of two standing waves of frequencies u> and 2uj: 

u(x, t, t) = p(x, t) suit + q(x, t) sin 2r (7.26) 

£(x, t, t) = —p(x, t) cos r — — q(x, t) cos 2r (7.27) 

[u,$] = i [p, q] (2 cos t sin 2t — cos 2r sin r) (7.28) 

Hence (|4T20|) yields 

^0 = i([«i]> = 0, Vi = l{[[u,tU}) = ll\p,q},p] (7.29) 

These expressions produce infinitely many examples of flows with 0(l)-drift in a super- 
critical asymptotic family (|5.15|) . 

7.6. A plane travelling wave of a general shape 
The velocity field (|7.1[) for a plane wave in three-dimensional space is 

u = Af'(k -x-r) (7.30) 

where A and fe are two constant vectors, / is a T-function of a scalar variable, its 
r-periodicity automatically leads to periodicity w.r.t. fc • x, primes stand for ordinary 
derivatives. Calculations yield 

Vo = A(A-k)(f' 2 ), Vx = A(A-k) 2 (f' 3 ) (7.31) 
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It shows that the only case of a zero drift corresponds to a transversal wave A_Lfc (an 
incompressible fluid) and the maximal drift takes place for a longitudinal wave A\\k. A 
pseudo-diffusion matrix is 

X lk =0 for (Mk)=AA k (f 2 ), (7.32) 
Next step is to consider the superposition of two travelling plane waves (|7.30D : 

u = Af(k ■x-t) + Bg'(l ■ x - r) = -Af T (k ■ x - r) - Bg T (l ■ x - r) 

\=-Aj T {k-x-T)-Bg T {l-x-r) (7.33) 

with constant vectors A, B, k, I. Calculations yield 

Vo = (A ■ k)A{j' 2 ) + (B ■ k)B(g' 2 ) + ((A ■ l)B + (B ■ k)A) (f'g') (7.34) 

which exhibits the interference (third) term. If /' = sin(fc ■ x — t) and g' = sin(l ■ x — t) 
then (jT3"4]) gives 

2V = (A ■ k)A + (B ■ k)B + ((A ■ l)B + (B ■ k)A) cos(fc - l)x (7.35) 
while Vi is too cumbersome to be presented here. Also 

2(66) = AiA k + BiB k + {A l B k + A k B t ) cos(fc - l)x 
For two mutually perpendicular longitudinal waves (A|| fc)_L(_B||Z) one can obtain 

%k = ~\{AiB k + A k B t ) [(A ■ k) 2 - (B ■ I) 2 ] sin(fc - l)x (7.36) 
which for the plane flow A = (A, 0), B = (0, B) gives 

X lk = - 1 -AB[{A-k) 2 -(B-l) 2 } ft ^sm(k-l)x (7.37) 

Here we have the eigenvalues Xi = —Xi of opposite signs, they are also oscillating in space. 
Once again (as in Sect. 7. 3), we deal with the sign-indefinite case of pseudo-diffusion. This 
example can be linked to plane sound waves; the drift (|7.31|) can give an addition to an 
acoustic streaming (see e.g. Lighthill (1978b) I. 

7.7. The drift for the polynomial velocity i6.18\) 

The expression (|6.27|) for the 0(l)-drift within a super-critical family contains the term 
given in (|6. 15|) 

V 12 = ([v,?j]) (7.38) 

which is built by two mutually independent oscillating functions. Such a functional free- 
dom allows to construct broad functional classes of drifts. However, in this case one should 
keep in mind that the degeneration condition Vo = —w (|6.25|) imposes a restriction on 
the coefficients of (|6 . 18|) . For example, one can chose v such that Vq = —w = 0. In this 
case the function v in (|7.38[) is not arbitrary. Alternatively, one can consider Vq = —w 
as the definition of w. In this case the functions v and rj in (|7. 38[) are indeed mutually 
independent. 

7.8. The Bjorknes configuration of two pulsating point sources 

This example is aimed to clarify the global structure of drift motion in an oscillating 
flow related to |Cook (1882")| |Hicks (1879)] |Hicks (1890)] and to show that interesting 
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oscillating flow, which arc different from waves, do exist. An incompressible velocity 
from the class of flows (|7.1j) in Cartesian coordinates x is given by 



y = x — l. I = const 



(7.39) 



which represents a superposition of two oscillating point sources. Calculations yield a 
rotationally symmetric w.r.t. the Z-axis field: 



y o= r , 2( ,aM a l 2 (|y| a »-N 2 y), V -(xxy) = 

(\x\\v\) 



(7.40) 



Let us introduce cylindrical coordinates (p,cj),z) with an origin at x = 1/2. The related 
ODEs (UTS) are 



p = Mpz, = 0, 



M // 2 2 2 

z + P 

2 V 4 ' 



(7.41) 



! = |I| 



M = -36 



p 2 + z 2 _ p /4 

\xWy\* 



\x\ 2 =p 2 + (z-l/2) 2 , \y\ 2 =p 2 + (z + l/2) 2 

The qualitative dynamics of particles described by this system is: (i) two singularities 
at the points x = ±1/2; (ii) any point of the sphere p 2 + z 2 = I 2 /A represents an 
equilibrium; and (iii) the directions of particle motions inside and outside of this sphere 
are topologically opposite to each other. 



7.9. The complexity of particle dynamics for Vq 
Now we illustrate the complexity of particle dynamics (18.31) for the zero-order drift ve- 
locity V (EPOj) ( see 



Aref (1984)[ |Ottino (1989)] |Samelson fc Wiggins (2006)] ). Let an 
incompressible velocity (|7.1[) be 





f cosy \ 




p= 


o 


-( 









asm z 
b sin x + a cos z 
bcosx 



(7.42) 



where (x, y, z) are Cartesian coordinates, a, b arc constants. Either of this fields, taken 
separately, produces simple integrable dynamics of particles. Calculations yield 



V = 



— a sin y sin a; 
b sin z sin y - 
b cos z cos y H 



- 2b sin y cos z 
a cos x cos y 
2a sin x cos y 



(7.43) 



Straightforward computations for this steady flow exhibit chaotic dynamics of particles. 
In particular, positive Lyapunov exponents have been observed. This example shows, 
that the drift created by simple oscillatory field can produce complex dynamics. Since 



the averaged dynamics is chaotic then related results by Arnold (1964) Aref (1984) 
Ottino (1989)| |Samelson fc Wiggins (2006)| |Chierchia fc Gallavotti (1994)| can be ap- 
plied to it. This example brings up new questions: (i) what is the relationship between 
chaotic motions for the original dynamical system and the averaged one? (ii) can the os- 
cillatory part of a solution also cause chaotic dynamics? (iii) can a chaotic drift and the 
pseudo-diffusion of Sect. 4 complement each other? (iv) how a chaotic drift can be used in 
the theory of mixing? This example has been constructed and computed by A.B.Morgulis 
(private communications) for the use in this paper. 
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7.10. The appearance of pseudo- diffusion (PD) due to slow-time-modulations 

Let us consider a simple example that clarifies the meaning of pseudo-diffusion. 

An infinite fluid oscillates as a rigid body with the small displacement of any material 
particle x — x(t,r). The Eulcrian coordinate of a particle is 

x = x + x, (x)=0, x = ext, xi G O(l) (7.44) 

where (•) is the r-average (|2.17[) ; the small parameter e and the function x±(t, r) can 
be either chosen by us or taken from (|8.14[) - (|8.16[) . A drift for a rigid-body oscillations 
is absent V = 0, hence the mean coordinate sofa particle is not changing with time. 
For simplicity we also accept that the mean coordinate x coincides with the initial (La- 
grangian) coordinate X of a particle, hence x(0, 0) = 0. 

A distribution of a Lagrangian tracer a is given as a = A(X). Then 

a(x,t,T) = A(x - exi) (7-45) 

One can expand both sides of (|7.45[) as 

2 ^ , _ dA(x) e 2 d 2 A(x) 

a + eat + e a* + ■■■= A(x) - ex u ^ + -J XliXlk dx dx k + " ' ^ ^ 

where a n = a n {x, t, r) = a n (x, t) + a n (x, t, r). The average of this equation yields 

£ 2 d 2 A(x) 

a Q + eat + e 2 a 2 + . . . = A(x) + — — + ... (7.47) 

where we have changed Xt(t, r) to €(t,r) (|4.16p . which is valid for the given precision. 
The t-differentiation of (|7.47p gives 

a ot + ea u + e 2 a 2t + ... = i-(^ k ) t d + ... (7.48) 

where we have used A = ao (which follows from (|7.47[) ) . One can see that eqns. (|4.17l) - 
(f47T9]l taken for zero drift V = V + eVt + e 2 V 2 + . . . = coincide with (fT38]l . 

It means that in this particular case PD represents a diffusion-like term that appears 
in the averaged equations as a correction a 2 caused by the curvature of the function 
ao = ao(x). In order to make this explanation clear one can imagine a one-dimensional 
case of (|7.44|) with A = A(Z) in (|7.45j) for a single Lagrangian coordinate Z. Then the 
averaging of oscillations of the entire graph A = A(Z) in Z-direction evidently produces a 
second-order correction (proportional to the curvature Azz) to the distribution ao(z) = 
A(z) which occur in the absence of oscillations. Now one might assume that for general 
oscillating flows the nature of pseudo-diffusion is the same: 

Conjecture: Pseudo-diffusion in (|4.19|) . (|4.25[) is always caused by the oscillations of 
Lagrangian tracer a{X, t,r) with respect to fixed Eulerian coordinates x. 

According to this conjecture the presence of a drift produces only a logically natural 
change from d/dt (for flows (|7.45[) , (|7.48[) without a drift) to the 'material' derivatives 
(with the drift velocity V in (j47T7j)-(j47T9j)). Hence this conjecture looks reliable, it indi- 
cates that pseudo-diffusion represents a natural as well as necessary term in the averaged 
equations. 

Remark: The effects of pseudo-diffusion and diffusion (or anti-diffusion, or anisotropic 
diffusion-anti-diffusion) are described by the same equations, hence they are mathemati- 
cally equivalent to each other. However, pseudo-diffusion is used only with regular asymp- 
totic procedures. 
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7.11. The sub- critical family of oscillations: procedure with a = 1 

Let us consider an example of a sub-critical asymptotic family with a = 1 and f3 = 0. 
For this family [/ = L/T = const in (|2.13l) which might be seen as an advantage for some 
applications. Following (|3.10[) we accept that (|2.33|) is: 

u(x,t,r) = u(x,t,r), S = (7.49) 

Then ([234) yields: 

Da = wa T + a t + (S • V) a = (7.50) 

The small parameter ei = 1/oj (|2.9p allows to rewrite it as 

£»ia = a r + e(a t + tt- V)a = 0, Di = e£ = D /u (7.51) 

where the subscript in E\ has been dropped. We are looking for the solution of (|7.5ip 
in the form of regular series (|4.4[) with redifmed e. The substitution of (|4.4[) into (|7.51|) 
produces the equations of successive approximations 

a QT = (7.52) 

a nT = -(d t +u ■ V)a n _i, d t = d/dt, n = 1,2,3,... (7.53) 

The solving of ([732]) . ([733")) yields 

So = 0, (7.54) 

oi = -(€-V)a 0l (7.55) 

a 2 = -(?-V)ai -{(«-V)ai} T -aI tJ (7.56) 

aot = 0, (7.57) 

ait + (V • V)a = (7.58) 



a 2t + (Vo ■ V)oi + (Vi ■ V)a = ^ ( Xifc^ I , (7-59) 

vf = Vi + i([I,I t ]> + \(£divt) t , 2xT = (Mk)t, (7.60) 

where Vo and Vi are the same as in (|4.20|) . Equations (|7.57j) - (|7.59j) can be written as a 
single advection-pseudo-diffusion equation (valid with the error 0(e 3 )) 

fr+^-K^sl) (7 - 61) 

F = £ F„ + £ 2 V^ b , «K b = t 2 af (7.62) 
a = ao + eat + e 2 ci2 (7.63) 

Eqn. (|7.61|) shows that the averaged motion of a represents a pseudo-drift with velocity V 
and pseudo- diffusion with the matrix-coefficient k,^. The reason for introducing the term 
pseudo-drift is the same as we used for pseudo-diffusion in (|4.25[) : the drift-like terms in 
(I7.58l) . (|7.59p . ()7.6ip play a part of sources, known from previous approximations. 
Remarks to Sect. 7.11: 

1. A simplified (in comparison with (|4.23p ) expression for pseudo-diffusivity in (|7.59l) 
complies with the interpretation of pseudo-diffusivity of previous subsection. 

2. The system (|7.57[) - (|7.60[) produces mainly diverging solutions. Indeed, (|7.57p yields 
So = ao(x). Let clo(x) 7^ const, then (|7.58p gives a linear with t growth of ai for any t- 
independent Vo, for example for the Stokes drift (see Sect. 7. 3). If we choose ao = const, 
then a similar growth appears for Hi etc. A statement that any sub-critical asymptotic 
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procedure with at £ O(l) produces diverging solutions might be considered as a conjec- 
ture, but its analysis is beyond the scope of this paper. 

3. This equations of this subsection do not produce the transport equation (drift) in 
any approximation. 

4. The arguments of items 2 and 3 emphasize the significance of critical and super- 
critical asymptotic procedures (of Sects. 2-7). 

Remarks to all Sect. 7: 

1. From the examples of Sects. 7. 1-7.9 one can see that the expressions for the drift 
velocity V (|4.20[) , (|4.2ip . (|6.17|) can provide arbitrary functional form and magnitude not 
higher than 0(1). In order to make further progress one should prescribe some particular 
oscillating flows (|2.2j) . 

2. These examples also show that pseudo-diffusivity in (|4.25[) . (|4.26j) can appear in 
three qualitatively different cases: (i) all positive eigenvalues %s corresponds to ordinary 
diffusion; (ii) all negative \i corresponds to anti-diffusion, and (iii) the mixed signs cor- 
respond to more complex anisotropic evolution with diffusion in some directions and 
anti-diffusion in the others. The complexity of the problem (|4.25[) increases if one con- 
siders the dependence of Xi on x an d t. The general theory of related linear PDEs (e.g. 
([739]) ) has not been developed yet, see |Polyanin (2002)| [Myint-U k Debnath (1987) 



3. The number of our examples is naturally restricted, therefore Sect. 7 illustrates our 
studies of Sects. 3-6 only partially. 



8. Links to Other Theories 

8.1. TTAM-solution of characteristic equation 

We are not aware about any paper which calculates several approximations for a drift 
by solving the related ODE by two-timing method. Therefore we present the TTAM- 
versions of such calculations here. Let us return to the original dimensional equation 
(|2.1[) . (|2.2p with an oscillating velocity 

u* = u*(x*,t*,r) (8.1) 

where x* is an upgraded notation for original Eulcrian coordinates which is introduced 
here instead of x* (see Sect. 2) since Eulerian coordinates in Lagrangian description also 
represent the hat-functions (|2.15|) . The oscillating trajectories 

x* = x*(X*,t*,r) (8.2) 

(which represent the characteristics for the hyperbolic equation (|2.ip , (|2.4|) ) can be found 
by solving Cauchy's problem for an ODE 

^- = u*(x*,t*,r), x*\ t = = X* (8.3) 

where X* is the Lagrangian coordinate of a fluid particle and t* = automatically leads 
to r = 0. Using the same procedure as for (|2.10[) we can obtain the dimensionlcss form 
of (1831) 

^ ^ dx d d d 

x = x(t, t) : — — = ujoutx, t,T); t=s,r = ujs, — = — + uj— (8.5) 
ds ds ot or 

which represents the two-timing form of a dynamical system for the motion of particles, 
see |Arnold (1964)| |Aref (1984)[ |Ottino (1989)] |Samelson k Wiggins (2006)[ We accept 



26 V. A. Vladimirov 

(I2.33|) and introduce a test-solution for our inspection procedure 

x{t,r) = x (t) + — x x {t,T), a> 0; 6 = u fj -\ (3 < 1 (8.6) 

that is motivated similarly to p. II) (if one takes Xq ^ 0, then r ceases to be a fast 
variable). The substitution of (|8.6[) into (|8.5f) yields 

w ^ + 1) (*° + i s = ^ (*° + i 2 (8j) 

The decomposing of the RHS into Taylor's series with the retaining of two leading terms 
gives 

x ot + u x - a x lT + uj~ a xit = uj p u + cjP~ a (xi • V)S (8.8) 

where u = u(x(t, r), t, r). The bar- and tilde- parts of this equation are 

x ot + u- a x lt = w"- a ((5i • V)u) (8.9) 
uj l - a x lT + uj- a x lt = uj^u + u p ~ a (xi ■ V)u + u> p ~ a {(xi ■ V)u} (8.10) 

where we have used the brace notation (|2.28p . The leading terms in these equations are 

x ot = c^- a ((5i • V)u) (8.11) 

UJ^Xxr = UJ P U (8.12) 

Here the T-average (■) docs not carry any superscript, since the spatial independent 
variable is absent; at the same time u depends on an unknown function xo(t). The 
equation (|8.12|) gives a + ft = 1, after that (|8.11j) brings us back to the same notions 
of super-critical, critical, and sub-critical asymptotic families (|3.8|) - (|3 . 10|) . For a critical 
family a = j3 = 1/2 and e = ex m = which produces an asymptotic problem 

d d\ ^ ^ 1^1^ /n oN 

lu— + — ) x = ^uju, x = x Q H F a;H xi + ... (8.13) 

OT Ot J y/CJ U! 

or 

x T = —e 2 xt + su(x, t,r), x = x~o + exx + e 2 x~x + • • ■ (8-14) 

The successive approximations are 

x 0T = (8.15) 
xxt = u (8.16) 
x 2t + x ot = (xx ■ V)tt (8.17) 

S 3 r + Xxt = (X2 ■ V)tto + -zXxiXxk a_ 5= (8-18) 

2 dxoidxok 

X iT + X 2t = (X 3 ■ V)M + XxiX-Zk——— h -XxiXxkXlj „_ 7~ „_ (8.19) 

dxoidxok 6 axoiOXokOXoj 
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where uq = u(xo,t, r). The system (|8.15|) - (|8.19l) can be solved using the same method as 
(|4.5[) - ()4.7[) and (|8.56l) - (|8.59l) . Here we only present the solution for the averaged motion 

x t = Uo (8.20) 

X u = (X! • V)C7 + Vi 

x 2t = (xi ■ V)Ui + (x 2 ■ V)L/ + U 2 

U Q = ((x 1 -V)u), (8.21) 
U 1 = -(( Xl • V)(5 • V)Si) = i([[u,ii],xi]) (8.22) 

— m ^ / d 2 u(x)\ 1/ _ d 3 x(x) \ /oooN 

U 2 = ((* 3 • V)«) + ( x U X2 k —y- )--( XliXlkUlj —-±-L ) (8.23) 



which can be rewritten as 



x t = U + (I ■ V)U (8.24) 
x = x + ex i + e 2 Xi + . . . , T = x — x~o 



U = U + eU 1 +e 2 U 



2 



One can see that the system (|8 . 24[) describes two kinds of motion: in one every material 
point moves with velocity C7; in the other every small material arc SI is stretched by the 



same velocity field according to a standard law, see Batchelor (1967) 



.2. TTAM -approach to GLM-kinematics 



The most general formula for a drift was obtained in GLM-theory by Andrews & Mclntyre (1978) 
which is referred hereafter as AM. The direct comparison of the results of Sects. 2-7 with 
those of AM is not feasible due to the following conceptual differences: (i) GLM-theory 
in AM is not fully adapted to the two-timing method; (ii) AM uses one small parameter 
(the amplitude of 2;) while we operate with two small parameters; (iii) AM employs a 
different averaging operation; and (iv) AM expresses a drift in terms of generalized La- 
grangian displacements x (|8.37p . not a given velocity field (as in our case). We develop 
here the TTvlM-version of GLM-kinematics with the aim to compare it with the results 
of Sects. 2-7. 

At the beginning of this subsection we use dimensional variables, but for brevity we 
suppress the asterisks; the use of dimensionless variables (starting from (|8.48p ) will be 
additionally notified. The core of the whole GLM-theory is a transformation of (|8.3[) into 
a PDE. In order to perform this transformation, we switch to mutually independent t, t 
(see Comment C to (|2.16|) ) and introduce a r-averaged trajectory (|8.2[) as 

x = x(X,t) = (x) x (8.25) 

where (-} x emphasises that the r-average (|2.17|) is performed for fixed X. From (|8.3j) we 
get x(X,0) = X. Equations (|8.3|) . (|8.25[) give us three pairs of main kinematical items: 

(i) two kinds of trajectories for a selected particle 

x = x(X, t, t) - the family of exact trajectories, (8.26) 
x = x(X, t) - the averaged trajectory, (8.27) 
with x(X,0,0) = x(X,0) = X, (8.28) 

(ii) two kinds of velocities 

u = dx/ds\x - the family of exact velocities of a particle, (8.29) 
V = dx/ds\x - the averaged velocity of a particle 
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(iii) two Eulerian expressions for the same material derivative 

d/ds\ x = d/ds\ s + u ■ V = d/dsfe + V ■ V (8.30) 

Relations (|8.26l) - (|8.30l) show that for Lagrangian description (X,t) one can use two differ- 
ent Eulerian descriptions (x, s) and (x, s) where x and V arc called the GLM- coordinate 
and GLM-velocity. One also should keep in mind that in (|8.29[) . (|8.30[) according to the 
chain rule 

d/ds\ x = d/dt\ x . T + u d/dr\ x , t (8.31) 

and similar equalities for fixed x or x. The existence of one-to-one mapping between X 
and x represents a key postulate of classical fluid dynamics; it guarantees the existence 
of the unique inverse to (|8.26|) function 

x = x(X, t, r) «• X = X(x, t, r); < J < oo, J = ' y\ (g 32) 

with a Jacobian Jq. After that the function 

x = x(x,t, t) (8.33) 

can be obtained by the completely legal substitution of X = X(x,t,r) (|8 . 32[) into 
x = x(X,t) dHHJ). 

The additional assumption (specific for GLM) is an assumption of invertibility of (|8.27l) 

x = x{X 1 t)^X =X(x,t); 0< Jr < oo, J x = (8.34) 

The coordinate S is defined by the averaging operation (|8.25|) . hence it does not represent 
any physical motion; (|8.34l) represents an assumption which can be valid only under some 
additional restrictions. For the GLM-theory this invertibility is absolutely required; let 
us also accept it and show how it leads to GLM-kinematics. 

The substitution of the inverse function X = X(x, t) (|8.34p into x — x(X, t, r) (|8.26|) 
produces the key (for the GLM-theory) relation between x and x 

x = x(x, t, r); < J 2 < oo, J 2 = d( f 1 ^ 2 ^ ( 8 _ 35 ) 

(#1,2:2,2:3) 

that is inverse to ()8.33j) . The functions (|8.33j) . (|8 . 35[) give us one-to-one mapping between 
two systems of Eulerian coordinates x and x; since J2 = Jo-/] -1 , the conditions for the 
Jacobians in (|8.34[) and (|8.35|) are equivalent to each other, provided (|8.32[) is valid. 
An important property of the averaging operation (|8.25p 

(•r = X (8-36) 

follows from the fact that r is not involved into transformations (|8.34j) . (|8.25D between 
X into x. In other words, the r-averaging operation commutes with the changing of 
variables (JSHJ). Taking (-) w of we get 

x(x,t,r) = x + x(x,t,r) with (x) x = x, (x) x = (8.37) 

where (x) x = x is valid by virtue of (|8.25[) and (|8 . 36[) . Simultaneously (|8 . 3T[) represents 
the definition of x which automatically possesses zero r-average; in AM x is called 
generalized Lagrangian displacement. 

Applying two Eulerian forms of material derivative (|8.30[) to both sides of (|8.37[) we 
obtain 

(d/ds\ s + u- V)£ =u(x,t,T) = (d/ds\ w + v- V){x + x(x,t,r)) (8.38) 
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which can be rewritten as 

d d — \ 

— +Lu—+V-Vj(x + x(x 1 t,T))=u(x + x(x,t,T),t,T) or (8.39) 
d d — \ 

— + w— + V ■ V x(x, t, t) = u(x + x(x, t, r), t, r) - v(x, t) (8.40) 

Ot OT J 

where (|8.39[) represents the equation for characteristics (|8.3[) written as a PDE in vari- 
ables x, t, t. We call either (|8.39j) or (|8.40j) the two-timing form of Andrews- Mclntyre- 
Kinematics-Equation (AMKE); it contains unknown functions x(x, t, r) and Ts(x, t). The 
description of fluid kinematics by (|8 . 39[) is rather unusual: AMKE describes an effective 
medium which moves with the averaged velocity V that includes both an advective ve- 
locity and a drift. The exact motion of a material particle in this medium is described 
by two different fields: the averaged motion is described by GLM-velocity V(x, t) and 
oscillatory displacement (from the GLM-position x) is given by x(x, t,r). 

The first step in solving (|8.40[) is straightforward: its bar-part (-) x gives us GLM- 
velocity V expressed in the terms of an original velocity 

T7(x, t) = (u(x + x(x, t, r), t, t))*, (8.41) 

which represents the main kinematical formula by AM. The average operation in the RHS 
of (|8.41[) is called GLM- averaging. GLM-averaging of a function a(x,t,T) is denoted by 
AM as 

a L = a L (x, t) = (a(x + x(x, t, r), t, r)) ar (8.42) 
The advection of a Lagrangian marker is described by the equation (|2.1[) 

{ d/ds\ B + u(x, t, t) ■ v) a(x, t,r)=0 (8.43) 



Due to (|8.30l) . we replace a material derivative in (|8.43p with d/dt\x + v ■ V and use 

(d/dsy + V ■ V)a(x + x(x, t, r), t, r) = (8.44) 

Then (-)^ yields 

(d/dty+u- V)d L = (8.45) 

This remarkable equation says that the GLM-average a L (|8.42p in the GLM-coordinates 
x ([8725]) . (|8727| is purely advected with the GLM-velocity V ([OTjl . The drift by AM is 

U(x, t) = V(x, t) - (u(x, t, t)) w = (u(x + x{x, t, r),t, r)f - (u(x, t, T )) w (8.46) 

For a purely oscillating field u = u the average (u) x = 0, hence 

U = u(x, t) = (u(x + x(x, t, t), t, r)f (8.47) 

To follow the route used by AM we can take a; of a small amplitude and decompose (|8.46[) 
and (|8.47p into Taylor's series; however in the two-timing case this problem contains two 
small parameters and is richer in asymptotic procedures. In order to describe them we 
write the dimensionless form of AMKE (|8.39|) which is obtained by the same procedure 
as ([2~T0| 

d d — \ 

— + lu— + v ■ V )(x + x(x, t, r)) = loSu(x + x(x, t, r), t, r) (8.48) 

Ot OT J 

where the given velocity field is taken as a purely oscillatory one. Starting from (|8.48p 
and further in this section we use only dimensionless variables and functions, keeping the 
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same notations without asterisks. Using the same inspection procedure as in Sect. 3 we 
accept (|2.33|) and introduce a test-solution (similar to (|3.ip ) 

x = x + x = x H — y(x, t,r), a = const > (8.49) 

where 16 BnO(l), y 6 TflO(l) and the amplitude of an oscillating part is given by 
the small parameter e a = l/uj a (|2.9p . Hence we obtain 

-| + u^. + V • + ^ V(x, t, T )j = oj^u + ^ V(x, t, t)\ (8.50) 

Decomposing the right-hand side of (|8 . 50[) into Taylor's series and omitting the terms of 
order 0(l/ui 2a ) and above, we have 

V + u l - a y T + = ujPu(x, t, t) + Lo?- a (y ■ V)u(x, t, r) (8.51) 

The bar-part of this equation is 

u = u f) - a {{y- V)u) (8.52) 

while the difference between (|8.51[) and (|8 . 52[) produces the tilde-part of the equation. 
To make a meaningful equation out of the tilde-part (cf. with Q3.4p ) it is necessary to 
accept that dominating terms are of the same order: 

a + /3 = 1 and y T = u(x, t, r) (8.53) 

Equations (|8.52p and (|8.53l) brings us back to the same notions of super-critical, critical, 
and sub-critical asymptotic families (|3.8[) - (|3.10[) where a critical family is characterised 
by a = (3 (a classical AM case corresponds to a sub-critical family with a = 1 and 
/3 = 0). Then AMKE (|8.39j) for the critical asymptotic family a = ft = 1/2 takes form 

d d — —\ 

w— + — + [/• V (x + x(x, t,r)) = s/uu(x + x(x, t, r), t, r) (8.54) 
or at J 

xCx, t, t) = —=xi(x, t, t) H — x 2 (x, t,r) + ... 
U(x, t) = U {x,t) + -r=Ui (x, t) + -U 2 (x,t) + ... 



The use of e = e 1 / 2 = 1/a/w (|2.9|) transforms (|8.54|) into 

'■(Uo + eUi + e 2 U 2 + ...)■ VJ (x + ex! + e 2 x 2 + . . .) 



d_ 2 d_ 

= eu{x + exi + e 2 x 2 + ...,t, r) (8.55) 

The first four approximations of this equation are: 

xi r = u(x, t, t) (8.56) 

x 2t + U q = (x! ■ W)u(x, t, t) (8.57) 

\ d 2 u(x t t) 

x 3t + x u + U x + (U ■ V)Si = (x 2 ■ V)u(x, t, t) + -x u x lk gl^J (8.58) 

xat + x 2t + U 2 + (Ui • V)xi + (C/ • V)x 2 = (8.59) 
,~ m-v- . x , ~ ~ d 2 u(x,t,T) 1„ „ „ d 3 u(x,t,r) 



oxiOXk o oxiOXkOx 
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These equations can be solved in the same way as the equations in Sect. 4- The solutions 
are 

xi = u T (8.60) 

x 2 = {(2i • V)S} T (8.61) 

— 1 f d 2 u(x) 1 T 

x 3 = x\ t + (U ■ V)x\ + {{x 2 ■ V)«} T + - j x u S lk g _ ^_ J | (8.62) 

with the same U ,Ui, and U 2 as in fOIjl - pT^ . 
Remarks and discussion: 

1. In (|8.48l) we are looking for the dimcnsionlcss drift of order one V = 0(1) that (as 
we know from Sects. 4 and 5) is true for critical and super-critical oscillations. The more 
systematic (but more cumbersome) approach is to replace V in (|8.48|) by ujSv (which is 
required by (|8.47[l ) and to derive rigorously that the first non- vanishing term in uj5V is 
0(1). 

2. The super-critical versions of the problems considered in Sects. 8. 1,8. 2 can be solved 
similarly to the problems in Sects. 5, 6. 



3. The expression f8T20]) . (|OT]l was obtained by |Yudovich (2006) 

4. The theories of both Sect. 8.1 and Sect. 8. 2 operate with the average taken for a 
fixed Lagrangian coordinate X . Indeed, it is the only possibility for the problem (|8.5j) : 
this problem describes dynamics of a single particle. As to AMKE (|8.39p . it explicitly 
uses the average for fixed x which is the same as for fixed X, see (|8.36[) . 

5. One can observe a strong resemblance between the system (|8.15[) - (|8.19p and the 
system (|536>(|g39")) . In (|5T5|) - (|5rg)) terms with U ,Ui,... in LHS are absent and 
Taylor's series in RHS contain x = x + x (instead of a; in (|8.56[) - (|8.59|) V These differences 
are due to the fact that x is an independent variable in (|8.56[) - (|8.59[) but in (|8.15j) - (|8.19j) 
it play a part of an unknown function. However one can show that these two systems 
are mathematically equivalent to each other, and the drift velocity U in Sect. 8.1 and 
in Sect. 8. 2 is the same. These properties are natural, since AMKE (|8.39|) represents a 
transformed ODE for characteristics (|8~3"1) (sec (|8T25j) - (J8T39J) ) . 

6. Two sets of formulae for drift velocities (|4TTC|) . (|4T2U)) and (f8TT5U|) . (|8~^T|) . (jg^j) look 
identical if one makes a correspondence 

u(x,t,r) -s-> u(x,t,r), £(x, t, t) O x(x, i, r), (8.63) 
V (x,t) o U (x,t), Vi(x,t) «-» t7i(sc,i) 

hence one may conclude that our formulae for the zeroth and first approximations of 
a drift are the same as AM. However this resemblance is misleading, at least partially. 
The crucial point is the use of the GLM-coordinates x (|8.27[) in GiM-kinematics vs. the 
original Eulerian coordinates x = x (|8.26[) in Sects. 2-7. It leads to the different definitions 
of the averaging operations with fixed x vs. fixed x. For example, the expressions x\ = u T 
(|8.60[) and £ = u T (|4.16[) look identical, however they are different since the former is 
expressed as a function of x and the latter as a function of a; = x: these variables are 
constant during the T-intcgrations. The implications of such a difference lay beyond 
the scope of this paper. At this stage one can see that: (i) the similarity between the 
expressions in (|8.63|) indeed corresponds to the close results for the zeroth and first 
approximations (due to the smallness of the difference x = x — x); (ii) the expressions 
for U2 and V 2 do not exhibit the similarity shown in (|8.63|) . 

7. One can see that, in order to express (|8.45j) in the spirit of Sects. 2-7, a L has to be 
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decomposed into Taylor's series for small x. The resulted equations is 

{d/dt\ w + Ts-V)R=0, (8.64) 

-pr , ..=• /_ 9a(x, i, r) \ x 1 / d 2 a(x,t,T) 

R= (a(x,t,T)) x + {x t v a l 7 ) + -(0*0* 



95;, / 2 \ dxidxk 

where x(x, t, r), X/(aF, t), and a = a + a must be expressed as in (|8.54[) and (|4.4j) . Then 
the terms of the type (xida/dxi) should be calculated with the use of the full equation 
(|8.44[) . Such calculations will produce Ricmann's invariant R which is constant along the 
averaged characteristic curves (or R is transported with the drift velocity). 

8. It is apparent that the establishing of one-to-one correspondence between (|4.25[) and 
(|8.64[) is not feasible. Indeed, (|4.25[) has a 'dissipative' pseudo-diffusion term which can 
not be incorporated in the transport equation (|8.64[) . This fact should be expected from 
Sect. 7. 10, where we have shown that pseudo-diffusion appears due to oscillations of La- 
grangian coordinates with respect to fixed Eulerian coordinates. Hence, pseudo-diffusion 
(which is a product of Eulerian averaging) is not compatible with (|8.64j) (which is a 
product of Lagrangian averaging, see the item 4 above) . This incompatibility complies 
with a general understanding that different averaging operations preserve and loose orig- 
inal information differently. One can conclude that the averaged fields in TT.AM-form of 
GLM-theory and in the theory of Sects. 2-7 contain complementary (or just additional 
to each other) information. It opens the opportunity to look for the advantages given by 
the presence of these two descriptions. 

8.3. Links to other classical theories 
The classical drift velocity. The classical expression for a drift velocity is given by the for- 



mulae (26) of |Longuet-Higgins (1953)] or by (5.13.21) of |Batchelor (1967)| Our V lEOO) 
will be proportional to this expression if it is rewrittenjvith the use of T-periodicity 
and integration by parts. In fact, the conjugated form ((£ • V)5>) of Vq = —((v ■ V) £) 
(|A 17[) coincides with the classical expression. Consequently Vq gives the correct expres- 



sion for the drift in a potential travelling wave (Stokes (1847)). At the same time, there 
are four main differences between our results and that of the quoted classical papers: (i) 
we describe all involved variables and fields precisely, within the two-timing framework; 
(ii) we employ two small parameters instead of one; (iii) we obtain different magnitudes 
for a drift while in the classical paper the leading term is quadratic in small amplitude 
(which is also obtainable within our approach); (iv) we obtain two more formulae for 
drift velocities as well as find pseudo-diffusion. 

The Krylov-Bogoliubov method. The Krylov-Bogoliubov averaging method (KB AM), 
see |Bogoliubov fc Mitropolskii (1961)||Krylov fc Bogoliubov (1947)[|Sanders fc Verhulst (1985)[ 



is aimed to solve problems for special classes of ODEs with oscillating coefficients. This 
method can be straightforwardly used for the calculating of the averaged equations of 
characteristic curves (|8.5p . We have exploited this option: the required calculations are 
rather cumbersome, therefore we formulate here the results only. The first iT.£L4M-term 
for a drift velocity coincides with Uq (|8.20j) . while the obtaining of the next two terms 
requires substantially more KBAM analytical calculations than TTAM ones. Here we 
make two comments: (1) KBAM as well as TTAM calculations of Sect. 8.1 give the av- 
eraged equations for characteristic curves. Within these two methods applied to an ODE 
the problem of finding the averaged Riemann invariants cannot be addressed. (2) The 
areas of applicability of TTAM are much broader than KBAM: for example TTAM can 
operate with PDEs and incorporate molecular diffusivity to the governing equations (see 
Sect. 9). 
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The homogenization theory. TTAM has the same methodological roots as the method 



of homogcnization (see Bensoussan, Lions and Papanicolaou (1978) Berdichevsky, Jikov, and Papanikolaou (1999) I, 



which represents a version of the two-scale method. The long- wave 'homogenization' solu- 



tions for the transport of scalar admixture obtained by McLaughlin, Papanicolaou & Pironneau (1985) 



Vergassola & Avellaneda (1997) Frisch (1995) show the appearance of a 'turbulent' dif- 



fusion matrix, which is always positive-definite. In order to establish further links we 
consider the case when the given velocity does not contain high frequency oscillations 
and can be expressed as Ui(x,t/u>,t) (or just u = ui(x,t)) with the 27r-periodic de- 
pendence on t. Let us show that in this case the problem can be transformed to the 
considered in Sects. 2-7. 

Problem A: The dimcnsionless form of the transport equation (|2.10l) . (|2.34l) for a purely 
oscillating velocity 

»W5.v)a-o, + 4 (8.65) 

where a = a(x,t,r), u = u(x,t,T), f3 = const, the mutually dependent time-variables 
are t = s, r = ujs. 

Problem B: A similar equation for different time-scales is 

|.+ W *u 1 .vU = > °=l(° +u °) ( 8 . 66) 

os J OS UJ \Otl OTl J 

where ai = ci\{x, ti, n), tti = u(x,ti,Ti), (3\ = const, the mutually dependent time- 
variables are t\ = s/oj, ti = s. 

One can see that the equations (|8.65[) and (|8.66p arc mathematically identical to each 
other if one introduces the link l Problem A^Problem B' as 

i-H-ii, r-H-n, /3o/3i + l, a t , a T G O(l) a tl , a T1 € O(l) (8.67) 

After such replacements any solution of (|8.65[) simultaneously produces a solution of 
(|8.66[) . Hence this rescaling procedure delivers a counterpart a(x,ti,Ti) of a(x,t,T): in 
the Problem B (|8.66j) the imposed oscillations have the frequency 0(1) and the slow- 
time-scale 0(1/ uj). The link (|8.67[) allows us to add one more solution for each solution 
in Sects. 2-7. The main motivation behind the rescaling (|8.67p is: the only solutions con- 



sidered in the homogenisation theory, see Bensoussan, Lions and Papanicolaou (1978) 
Berdichevsky, Jikov, and Papanikolaou (1999) are of the type (|8.66|) . 



The presence of molecular diffusivity. Let fi* and /i be the dimensional and dimension- 



less coefficients of molecular diffusion. Eqn. (|2.ip describes the advection-diffusion of a 
scalar admixture in an incompressible fluid after adding the term /z*V* 2 a to the RHS. 
All results of Sects. 2-6 can be straightforwardly generalized to the flows with fi = 0(1). 
Such generalization is achieved by replacing operators d/ds* — > d/ds* — /i*V* 2 in (|2.ip 
and d/dt — > d/dt — /iV 2 in (|2.4|) and in all subsequent formulae The limit /j — > is a 
regular one for our strictly regular asymptotic procedures (|4.4p , so all the problems with 
/i have the cases studied in Sects. 2-6 with /i = as their limits. 

Magnetohydrodynamics (MHD). One can also develop the same asymptotic procedures 
and to derive the averaged equations similar to (|4.17p ~ (|4.22p for a vectorial passive admix- 
ture such as magnetic field h* ' (x* , s*) in the kinematic MHD-dynamo problem for a given 



oscillating velocity field, see Moffatt (1978) In this case one can start with equations 



dh*/ds* + [h*,u*} =0, divft,*=0, divw*=0 (8.68) 



which replace (|2.ip . Taking the same velocity field (|4.ip and repeating the same steps as 
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in Sect. 4 lead to the averaged equation (valid with the error 0(e 3 )) 

_ __ 8 / dh\ — 

h t + {V -V)h-{h-V)V = — \K lk — y divfc = (8.69) 



h = ho + ehx + e t h 2 , V = V + eV\ + rVj, n ik = e% k 

where all coefficients are the same as in (|4.17j) - (|4.22j) . For this problem we have per- 
formed the detailed calculations of Vo and Vi, while the expressions for V2 and Kik 
represent reliable conjectures. In these calculations one should essentially use (|2.24|) and 
(|2.32p . Notice that eqn. (|8.69[) describes evolution of h(x,t) for an arbitrary spatial 



scale L, while the homogenisation approach to the same problem (see e.g. Frisch (1995) 



Zheligovsky (2009) ) operates with long-wave solutions. However, the rescaling of spatial 



variables (similar to one for time variables in (|8.65j) . (|8.66[) ) is applicable. 



9. Discussion and Conclusions 

1. The motivation behind this paper is to study in full a variety of asymptotic solutions 
and procedures for the advection of a scalar or vector field in an oscillating flow. The 
main result of the paper is the general understanding of the differences between different 
classes of asymptotic solutions. 

2. The author is not aware about any research paper entirely devoted to the general 
analysis of motions of a scalar (or vector) admixture in high-frequency oscillating flows. 
At the same time this topic is exploited in many papers in the solutions of particular 
applied problems (e.g. Magar & Pedley (2005) I. We hope, that our general study can 
underpin further applied studies. 

3. Several well-known papers are devoted to the motion of a scalar admixture in spa- 



tially oscillating flows, e.g. McLaughlin, Papanicolaou & Pironneau (1985) Vergassola & Avellaneda (1997) 



and a review in Frisch (1995) These authors have considered only one asymptotic family 
of solutions (following to Bensoussan, Lions and Papanicolaou (1978)). Therefore it can 
be useful to study both high-frequency and short-wave problems at a greater level of 
generality and then to consider different asymptotic procedures for flows that oscillate 
both in space and time. 

4. In our paper we are dealing with the systematic calculations of solutions. A few 
restrictions have been accepted at the very beginning, the rest of the paper does not con- 
tain any physical input. Hence, our results require physical interpretations, explanations, 
and applications. 

5. The asymptotic validity of our results is provided by the fact that we build only 
the solutions for high frequency uj. In our approach oj is not linked to any dynamical 
equations. Instead, we require that dimensional frequency lo* is high in comparison with 
the inverted slow-time-scalc 1/T (|2.7j) of an arbitrary prescribed velocity field (|2.2p . If 
velocity field (|2.2[) is purely oscillatory, then the choice of T is not unique, which is 
explained in (|237)) , flUM) . 

6. All results of this paper have been obtained for the r-periodic functions from the 
class H (|2.15j) . which is closed with respect to all used operations. One can try to con- 
sider more complex classes of quasi-periodic, non-periodic, or chaotic solutions. The 
generalization of our results to quasi-periodic solutions (containing several r-periods) 
looks rather straightforward, provided one can deal with resonances. At the same time, 
chaotic oscillations (containing low u in their spectrum) can be hardly treated by our 
method since u enters the denominators of asymptotic series (the apparent requirement 
for the applicability of our method is separating the frequency spectrum from zero). 
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Some generalizations of our results can be achieved after replacing the integration over 
to < t < to + 2ir by integration over — oo < r < +00 in the definition of average (|2.17[) 
(as it has been widely accepted for the spatial average in the homogcnization theory, see 



Berdichevsky, Jikov, and Papanikolaou (1999) I 



7. In principle, TTAM allows to produce the approximate asymptotic solutions with 
as small an error (RHS-residual) as needed. However, the next logical step is more chal- 
lenging: one has to prove that a solution with a small RHS-residual is close to the exact 
one. This problem is equivalent to the presence of a small additional force. Such proofs 



had been performed by Simonenko (1972) Levenshtam (1996) for vibrational convection. 



Similar justifications of TTAM for other oscillating flows are not available yet. 

8. Different asymptotic procedures £1/2, £1/3, £1/4, etc. correspond to different asymp- 
totic paths on the plane of two scaling parameters (|2.12[) . Two additional asymptotic 
paths represent two successive limits: (i) first 5 — > and then l/u> — > or (ii) first 
1/uj — > and then 5 — > (these paths are not considered in this paper but they are 
worth studying). The relevance of different available paths to particular physical situa- 
tions represents a major problem for further studies. At this stage we can mention only 
that different asymptotic paths correspond to different relations between physical pa- 
rameters. For example, 5 = uj^ 1 ^ 2 (|4.ip and 6 = uj- 1 ^ (52]) in their dimensional form 

give 

U = LT- 1 ' 2 ^ 1 ' 2 and U = LT^u* 2 / 3 (9.1) 

correspondingly. For the asymptotic procedures with large oj* (|2. 13[) it means that the 
imposed oscillatory velocity (|2.2j) is higher in the latter (super-critical) case. 

9. To be able to compare different flows physically one might introduce a drift- efficiency 
of an oscillatory flow as the ratio 

E a = (amplitude of drift velocity) / (amplitude of velocity oscillations) 

with a defined in (|2.9p . The dimensionlcss drifts in both cases (|9.1[) arc 0(1) and in 
dimensional form they are of order L/T. Then one might conclude that the efficiencies 
of (|9.1[) are E 1 / 2 = (cjT) -1 / 2 and E 1 / 3 = (u*T)~ 2 / 3 correspondingly, which indicates 
that the critical family solution (|4.ip produces stronger drift than the super-critical one 
(|5.2p . However, this comparison is not fair since the solution (|5.2p is valid only when the 
leading term in (|4.ip vanishes (|5.ip . Taking the degeneration into account one obtains 
the efficiencies of (|9.ip as 

E( /2 = (w*T) _1 and E 1/3 = (w*T)~ 2/3 

where Ef, 2 represents the drift-efficiency of the degenerated critical solution (|4.ip . (|5.ip . 
Now, one can conclude that a super-critical solution is more drift-efficient than a critical 
one. 

10. The appearance of small pseudo-diffusion (PD) in the high-frequency asymptotic 
problem (|4.2p . (|4.17p - (|4.27p is an accurate and qualitatively new result. One can make 
two conjectures: (i) a solution can slowly self-concentrate (due to small negative pseudo- 
diffusivity) or it can undergo unusual anisotropic evolution (due to sign-indefinite pscudo- 
diffusivity) ; and (ii) the maximum principle can be violated (since the original equation 
(|2.ip expresses the conservation of a in each fluid particle, hence the values of sup a and 
inf a do not change with time). One can argue that the conjecture (ii) is not valid due to 
the explanations given in Sect. 7. 10, while (i) requires additional studies. 

11. It is worth introducing the finite molecular diffusivity fi — O(l) (see Sect. 8) 
that can improve the convergence of all used asymptotic procedures. The introducing 
of asymptotically small or large diffusivity (e.g. /1 = 0(1/lu) or /1 = O(oj)) will generally 
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lead to different asymptotic theories corresponding to three independent small parame- 
ters from the extended list J22 



12. In Sect. 7. 9 we have shown that the drift, caused by a relatively simple oscilla- 
tory velocity, produces chaotic dynamics of particles. This result leads to numerous new 
questions (see the end of Sect. 7. 9) and deserves a serious elaboration. 

13. In Sect. 2. 3 we have demonstrated that an infinite and continuous range of slow- 
time scales 1/lo* < T < oo can be used in an important case where slow-time t is absent 
in the expression for velocity: u = u(x, r). We have formulated this result with the aim 
to clarify the existence of many scales (which often causes confusions) and to show that 
an infinite number of similar (to the considered in Sects. 3-7) solutions are available. At 
the same time, this result can lead to studying the motions with simultaneous presence 
of several different scales T, and to developing multi-scale (triple-scale, etc.) theories. 

14. It is worth completing the calculations of Ricmann's invariant and (|8.64[) and to 
compare the related solutions for a with that of (|4.25j) . The aim is to find the advantages 
given by two complementary averaged solutions of the same problem. 

15. The version of TTA M- theory for a vectorial passive admixture (with the av- 
eraged equations (|8.69[) ) is linked with the problem of kinematic M-fTD-dynamo (see 
Moffatt (1978) ) and can bring new results. It is apparent, that for the majority of drift 
velocities V(x,t) the stretching of material elements will produce linear growth \h\ ~ t. 
At the same time, there are known examples with an exponential stretching of mate- 
rial lines with time, say, in the flows near stagnation points or in chaotic flows. These 
examples will provide the exponential growth of \h\. 

16. Our theory brings up a new possibility: one can find an oscillatory velocity u 
which produces any required drift field (for example, a drift in the form of ABC-flow). 
The results of Sects. 4-7 (including (|7.7j) . (|7.38p ) indicate that the solution of this problem 
is not unique. 

17. TTy4M-method and results have potential applications in the studies of a broad 
variety of flows. The discovery of drift motions with the velocities Vi = 0(1) and 
V2 = 0(1) can lead to new applications. The other advantages of our method, that can be 
exploited in applications, are: (i) the mathematical generality of the imposed oscillatory 
velocity fields (|2.2[) . ()6.ip : (ii) the wide class of scalings (|2.14|) : and (iii) the most straight- 
forward Eulerian average. In particular, our method can be useful for applications to 



biologically motivated flows possessing complex geometry (e.g. Hydon & Pedley (1993) 



Magar & Pedley (2005) Leptos, Guasto, Gollub, Pesci, and Goldstein (2009) ), for micro 



hydrodynamics (e.g |Leal (1980)[|Kim fc Karriba (1991)[|Hinch fc Nitsche (1993)||Moffatt (1996)[); 
for flow mixing (e.ff. |Ottino (1989)[|Carlsson, Sen, and Lofdahl (2005)||Carlsson, Sen, and Lofdahl (2004)[ 
Leal (2007) I, for the flow of blood and flows in lungs, see Pedley (1980) and for the pipes 



with oscillating walls, see Carpenter & Pedley (2001) The spreading of ash from recently 
erupted Eyjaffjallajokull volcano has made the research on the general area of the trans- 
port of an admixture in a fluid (which can oscillate due to various reasons) even more 
important. 
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Appendix A. Calculations for a = 1/2 of Sect. 4 

The zero-order equation (|4.5|) is: 

a OT = (AI) 

The substitution of ao = ~a~o(x, t) + ao(x,t,r) into (|A lj) gives ao T = 0. Its TP-integration 
(|2.21| produces a unique (inside the T-class) solution ao = 0. At the same time (|A 1[) 
does not impose any restrictions on ao{x,t), which must be determined from the next 
approximations. Thus the results derivable from (|A 1[) are: 

ao(x, t, t) = 0, Vao = oq(x, t); =aa + ao = ao(x,t) (A 2) 

The substitution of the truncated solution into the governing equation (|4.3[) produces 
the RHS-residual of order s: 

T) 2 a [0] = Res[0] = e(u ■ V) a + e 2 a t = 0(e) (A 3) 

The first- order equation (|4.6|) is 

a iT = -(£t • V) a (A 4) 

The use of So =0 (|A 2|) and a\ T = reduces (|A 4j) to the equation ai T = — (u- V) ao. Its 
T-integration ()2.21j) gives the unique solution for ai 

Si = -(|-V)a (A 5) 

where £ = u T (|4.16j) . Hence, ai and a^l are 

ai = ai — (^ ■ V) ao, a' 1 ' = ao + eai ; Vao(a;, t) and a,\(x, t) (A 6) 

The substitution of a^ into f|4.3[) produces the RHS-residual of order e 2 : 

D 2 a [1] = Res[l] = e 2 ((u ■ V) a x + d t a^) = 0{e 2 ) (A 7) 



Formulae (|A2|) and (|A5j) give (jlTTTj) and (|4T2|) . 

The second- order equation f (|4.7p for n = 2) is 

a 2T = -(u ■ V) ai - a ot (A 8) 

The use of (|A 2|) and a 2r = transforms (|A 8|) to 

cl2t = ~{u ■ V)Si - (u ■ V) ai - aot (A 9) 

Its B-part is 

aot = -((«-V)oi) (A 10) 

where we have used (a2 T ) = 0, ((u ■ V) ai) = 0, and (aot) = ^ot- The substitution of (|A 5} 
into (|A 10|) produces the equation 

5ot = ((«-V)(i-V))5o (All) 

which represents a version of inspection equation (|3.13l) systematically derived. One may 
expect that the RHS of (|A lip contains both first and second spatial derivatives of ao, 
however all second derivatives vanish. In order to prove it we introduce a commutator 

K = [€, u] = (S ■ V)€ - (I • V)u, (A 12) 

(«■ V)(C- V) - (i- V)(«- V) =if V. (A 13) 
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The B-part of (jATg)) is 

((« ■ V)(£ • V)) = ((£ ■ V)(S • V)) + K ■ V (A 14) 

At the same time the integration by parts (|2.29|) gives 

((u-V)(€-V)> = -<(€-V)(tt.V)>, u = l T (A 15) 

Combining (|A 14|) and (|A15|) we obtain 

{{u-V){l-V)) = l -K-V (A 16) 

which reduces (|A lip to the advection equation (|4.17p with 

V Q ^-({u-V)i)=-\(\k,u]) = -\~K (A 17) 

which coincides with (|4.20[) . The T-part of (|A 9p appears after subtracting (|A 10|) from 
(|A 9|) and the use of notation (|2.28[> : 

2 2 r = -(«• V)oi -{(«• V)oi}. (A 18) 

Its T- integration (|2~21~j) with the use of (|A 5|) gives (|4~T3l> : 

S 2 = -(rV)ai + {(M-V)(i-V)} r a , VSi (A 19) 

Hence, a 2 and aft can be written as 

02 = 02 + 02, a' 2 ' = a + £ (a\ - (£ ■ V) a^j + s 2 a 2 , Vai,a 2 (A20) 

where clq and a 2 are given by f|4. 1 T[) , (|A 19|) . The substitution of into (|4.3p produces 
the RHS-rcsidual of order e 3 

£> 2 S l2] = Res[2] = e 3 ((« ■ V) a 2 + d t (a x + ea 2 )) = 0(e 3 ). (A 21) 

The third-order equation ( (|4.7|l for n = 3) is: 

a 3r = -(it ■ V)a 2 - au (A 22) 

Its B-part is 

o lt = -((«■ V)a 2 >. (A23) 

The substitution of (|A 19|) into (|A 231) . the use of it = £ T , and the integration by parts 
(|2~29)) yield 

ait = ((« • V)(€ • V))ai + ((I • V)(u • V)(€ • V))a (A 24) 

where ((it • V)(^ • V)) has been already simplified in (|A 16[) . The second term in the RHS 
of (|A 24p formally contains the third, second, and first spatial derivatives of ao; however 
all the third and second derivatives vanish. To prove it, first, we use f)2 .30[) : 

((I • V)(S • V)(£ • V)) = -((it • V)(i • V)(€ • V)) - • V)(I • V)(5 • V)) (A25) 

Then we use (|A 12[) . (|A 13[) to transform the sequence of operators (£ • V) and (it ■ V) in 
each term in the RHS of (|A 25|) into their sequence in the LHS. The result is 
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As the result (|A 24j) takes form (j4J8|) with V (|A 17|) and 

V t = -{{I ■ V)(u • V)0) = uU]> = (A27) 

which gives (|420|) . The T-part of (|A~22|) after its T-intcgration gives (|4T4"f 

a 3 = -(£-V)3 2 -{(w- V)a 2 } T -5L, SI = -(| T .V)ao (A 28) 

where a 2 is given by (|A 191) . 

Hence, S3 and can be written as 

a 3 = a 3 + a 3 , a m = a + e(ai - (£ • V) o ) + e 2 a 2 + e 3 S 3 , Va 2 ,a 3 (A 29) 

where a , ai, a 2 , and S3 are given by (|4.17|) . (|4.18p . (|A 19[) . and (|A 28|) . The substitution 
of a' 3 ' into (|4.3p produces the RHS- residual of order e 4 : 

Srfl = Res[3] = <3(e 4 ) (A 30) 

Explicit formulae for residuals (similar to (|A3[) . (|A 7|) . (|A21[) ) for (|A30|) . (gUg]) and for 
all other considered cases can be calculated straightforwardly; however for brevity we do 
not present them. 

The fourth-order equation ( (|4.7j) for n = 4) is: 

a 4r = -(« • V)a 3 - a 2t (A 31) 

Its B-part is 

3» = -<(u- V)o 3 ) (A 32) 

The substitution of a 3 (|A 28[) into (|A32|) . a 2 (|A 19[) into a 3 (|A28j) . the integration by 
parts (|2~29)1 . and the use of (|4T6]) . (|AT7]) . flX27| yield 

((5 • V)aa) = (Vo • V)a 2 + (Fi • V)3i + ((£ ■ V)(£ • V))(F • V)a - (A 33) 
-<(£ • V)(£ t • V))a + Xao, where V)(« ■ V){(5 ■ V)(£ ■ V)} r ) 

The shorthand operator X (as well as the operators 2), 21. 05, £, and 3" below) acts 
on ao- The RHS of (|A 33j) formally contains the fourth, third, second, and first spatial 
derivatives of ao; however all the fourth and third derivatives vanish. In order to prove it 
we first rewrite X as 

X= ((£• V)(«- V)2J r ) where 2) = (u ■ V)(£ ■ V) (A 34) 

The use of and ([ATI]) . (gTgJl transforms (|A~34l) to 

2X= -2T+S+ -((£• V)(£- V))(K- V) (A35) 

21= ((€-V)(e-V)(u-V)(e-V)), ®= ((^ T -V)(M-V)(i-V)) 
Let us now simplify 21 and 25. For 05 we use (12.301) 

B = ((k t • V)(« • V)(£ ■ V)) = —{(k ■ V)(€ ■ V)(f • V)> - ((lT • V)(4 ■ V)(S ■ V)) 

To change (£ • V)(u • V) into (2 • V)(£ • V) in the last term we use (|Xl2]l . (|A13|) that 
yields: 

® = -i((^.V)(€-V)(€-V))-i?g.V 1 £=[lT,i?] (A36) 
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The operator 21 is simplified by the version of (|2.30[) with four multipliers 

2T =<(€ ■ v)(€ • V)(S • V)(€ • V)) = -((a • V)(« • V)(€ • V)(£ • V)) - (A 37) 

-((€ • V)(S • V)« • V)« • V)) - ((I • V)(£ • V)« • V)(S • V)) 

The multiple use of commutator (|A 12|) . (|A 13[) allows us to transform the sequence of 
operators (£ • V) and (u ■ V) in each term in the RHS of (|A 37|) to the sequence in its 
LHS. The result is 

2T=-i((K.V)(i-V)(I-V)) + iK 77 .V, P = (A38) 
Now, (jA~35|) . (|A~36|) . and (|X38| yield 

X = i(F • V)((€ ■ V)(* ■ V)) + !<(€• V)(£ • V))(K . V) - i^+K 77 ) • V 
The substitution of this expression into (|A 33[) . (|A 32|) gives 

a 2t + (Fo • V)a 2 + (Vi • V)3i - -(« + K 77 ) • Va + 7O0 - = 0, (A 39) 

o 4 

c = (F • v)((I • V)(£ • v)) - ((€ • v)(« • V))(F • V), = ((i • v)& • v)) 

Additional transformations of the last two operators in (|A 39|) yield 

- * = ■ v - \({u> • v)c + (I • v)«0 • v - \ { u>i k + u' k l)^- k = 

= \(\£M-V-4r (x ik -£-) + kidivu' + u'div£) (A 40) 



u'^^-IFo.f], Xi fc = ^Ka+^) (A41) 

The substitution of ([A 40[) into ([A 39|) leads to the equation for U2 (|4. 19|) where the 
formula (|4.23[) for \ik 1S obtained from (|A 41[) by the use of definition u . 
The T-part of (|A 31[) after its T-integration gives (|4.15l) 

54 = -(C-V)a 3 -{(w-V)S 3 } r -^ t (A 42) 

where 03 is given by (|A 281) . 



Appendix B. Calculations for Super-critical Families 

B.l. Calculations for a = 1/3 of Sect. 5 

The equations of the zeroth and first approximations (|5.5p . ([5.6[) are the same as (|4.5p . (|4.6p 
hence they have the same solutions 

ao = ao(x,t), ai = 01- (£• V)o ; Va , ai G B n O(l) (Bl) 

An essential difference from the w 1//2 -procedure emerges for the second approximation 
(|5.7|) where instead of an advection equation for oq (|4.17[) we obtain a compatibility 
condition 

(y -V)a = (B2) 

which represents an identity due to (|5.1|) . The T-part of (|5.7p produces the same 02 as 
in (|XT9|) . Hence 

S 2 =a 2 +S 2 , a 2 = -(I-V)ai + {(M.V)(i-V)} T a , Va 2 6lfl O(l) (B3) 
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The equation for the third approximation (n = 3 in (|5.8[) ) jointly with (|5.1|) produces the 
B-part of the equation 

B(M = <(€-V)(tt.V)(€-V))Bb (B4) 
with the same triple-correlation as in |A26l) . (|X27| ; hence 

^+Fi-v)ab = 0, where F x = u],€]) (B 5) 

The T-part of (|5.8j) after its T-integration gives a simpler expression 

S 3 = -(i-V)a 2 -{(S-V)S 2 } T (B6) 

with a 2 (|B 3p . Hence, a 3 and aft can be written as 

03 = 03 + 03, a' 3 ' = ao + e(ai - (I • V) ao) + e 2 a2 + £ 3 a 3 (B 7) 

where ao, S2, and 03 are given by (|B 5[l . (|B 3p . and (|B 6|) . while a\, a 2 and S3 are to be 
found from further approximations. Hence the cj^-procedure shows that in the presence 
of degeneration (|5.1[) a drift velocity remains O(l) and is given by the expression V\ (|B 5p . 
The equation for the fourth approximation ( (|5.8[) with n = 4) is: 

7ii r = — (u ■ V) S3 — ait (B 8) 

After the use of the same steps as in the w^-procedure we have B-part 

a u + (Vi ■ V)3i + (V 2 ■ V)a = 0, (B 9) 

with the same notations as in (|B 5|) . ([5T7|) . ([535]) . The T-part of (|B~8)) after its T- 
intcgration gives 

a 4 = -(£-V)a 3 -{(u-V)a 3 y -a[ t (BIO) 

with a 3 ([B6]) . 



B.2. Calculations for a = 1/4 0/ Sect. 5 

The equation of the zeroth, first, and second approximations (|5.24p . (|5.25[) , and (|5.26p are 
the same as (|5.5p . (|5.6p , and (|5.7p so they have the same solutions (|B 1|) - (|B 3p 



a = a + a = Oo(a;, i), a\ = a\ — (| ■ V) ao, a 2 = a 2 + a 2 (Bll) 

An essential difference from the w 1 / 3 -procedure emerges for the third approximation 
(|5.26p where instead of an advection equation for ao (|B 51) we obtain a compatibility 
condition 

(Fi-V)a = (B12) 

which represents an identity due to (|5.20p . The T-part of (|5.27|) produces the same a 3 as 
in (|B 6p . The equation for the fourth approximation f (|5.28p with n = 4) is: 

a 4r = — (w • V) a 3 — a t (B 13) 

After performing the same steps as in the a; 1 ^ 3 -procedure we obtain B-part 

aot + (V 2 ■ V)a = 0, (B 14) 

with the same notations as in (|5.17p . (|5.18p . The T-part of (|B 13|) after its T- integration 
gives 

a 4 = -(£-V)a 3 -{(u-\7)a 3 } T (B 15) 
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where S3 is given in (|B 6j) . 
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